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A THEORETICAL STUDY OF THE DYNAMIC ERORERTIES 
OF HELICOPTER-BLADE SYSTEMS 
By H. Reissner and M. Morduchov 


SUMMARY 


The work herein presented on a theoretical study of the dynamic 
properties of lifting rotors covers: 

1. The derivations of the angles of attack of the inflow, of the 

"blade-posit ion variables — pitch, flapping, and lagging — and 
of the aerodynamic and inertia forces acting on hinged heli- 
copter (lifting-rotor) "blades in the hovering and in the 
traveling states 

2. The development and solution of the equilibrium conditions of 

the blade system in the hovering and in the traveling states 

3. The development of the frequency, stability, and damping proper- 

ties of the hinged, sufficiently rigid rotor blades during 
hovering and traveling 

The method of solution of small oscillations about a state of 
simultaneous rotation and traveling has in this paper been carried through 
for only small speed ratios. 

This method was applied to four cases of diverse constraint conditions 
between the three angles of pitch, flapping, and lagging. The results are 
significant in regard to restoring force and to modej frequency, phase, 
damping, and amplitude ratios. The number of modes and of independent 
amplitudes is, of course, equal to the degree of freedom. Each mode in 
hovering corresponds to one frequency, but in traveling each mode consists 
of three frequencies of fixed amplitude ratios and fixed phase, differences 
but with only one free amplitude. However, the amplitude ratios of the 
two additional frequencies to the original amplitude are, in all cases which 
have been computed numerically, smaller than the speed ratio. The results 
for the four cases treated show marked advantages obtainable by appropriate 
kinematic constraints 'between the three angles in regard to safety against 
resonance, to damping, and to automatic adjustment. 
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IHPRODOCTION 


The problems of blade— angle control, stability of motion, resonance, 
and flutter of helicopter-blade systems have not as yet been fully treated. 
This fact seems to be substantiated by the presence of disturbances in 
present-day rotor operations which have not been fully explained. No 
complete theory encompassing the effects of different methods of articula- 
tion and angle control appears to be known. These problems are closely 
interrelated through their dependence on the dynamic equations of blade 
motion. 

Previous publications have dealt mainly with the performance of the 
helicopter or with the stability of the equilibrium of steady flight of 
the complete helicopter system with very special assumptions in regard 
to the blade and hub connections. 

It is believed that the problem of smooth operation of a helicopter 
must be attacked in a more general way, and for this reason must be 
divided into at least two partsj namely, 

(1) The free oscillation of the blade system about the different steady 

states of flight 

(2) The forced oscillations of this system caused by the reaction of the 

fuselage, by irregularities of torque, by gusts, by transition to 
another state of flight, by flying in a curve, and so on. 

Problem (l) again falls into several parts. ■ The first part deals 
with the conclusions which can be drawn from the results of the theory 
of small oscillations about a steady state of motion applied to a system 
of sufficiently rigid blades hinged to a driving hub. The rigidity of 
blades is sufficient to give the real behavior of the blade system if the 
natural frequencies of a blade, treated as rigid, are small in comparison 
with the lowest natural frequency of elastic vibrations of a blade. 

The second part, not treated in this paper, would have to deal with 
the superposition and interference of elastic vibrations, or what is the 
same, of elastic waves of the blades on or with the rigid-blade oscilla- 
tions in those cases when the blades are appreciably flexible. The 
analysis of this phenomenon would require the integration of the equa- 
tions of deflection — and twist — vibrations of the blades, under the 
action of the local aerodynamic and inertia (including centrifugal) 
forces and under the effect of the boundary conditions at the hinge and 
at the tip. 

The flutter problem, meaning the determination of the critical 
velocity at which damping coefficients become negative, so that self- 
excited oscillations arise, haB not been covered in this paper for two 
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reasons. The first reason is that the critical velocity at the low average 

y 

reduced ratio — of the "blade will, in general, not "be reached, especially 
HR 

since the blade is stiffened by the restoring centrifugal force. The 
second reason is that all other sources of instability and resonance should 
be removed first before going into this difficult problem more deeply than 
previous authors, who have simply applied results of straight-moving air- 
foils. It is not improbable that sometimes unstable oscillations appearing 
in stationary flow, with no phase differences between bending and twisting, 
have been mistaken for flutter. 

This investigation was conducted at the Polytechnic Institute of Brooklyn 
under the sponsor shi p and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


% 




x,y,z 


*l* y l* z l 

ft 

v 

7 

* - fit 
t 

s 

d 

* ~ dt 

* wL. 

a* 

0 


a 


right-hand Cartesian coordinates, fixed to hub of blade 
(see fig. 1) 

right-hand Cartesian coordinates, fixed to blade system 
(see fig. 3) 

angular velocity of h.u"b 

velocity of flight (called traveling) in any direction, for 
example, forward, sideways, or backward 

angle between plane of rotation and velocity vector V 

angular position of blade center line (if ** 0 when x— axis 
coincides with projection of v on plane of rotation) 

time, seconds 

maximum thickness of blade cross section 


angle between plane of rotation and zero— lift line of chord 
(called pitch angle), positive from y to z 

angle of attack 



k 
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0 

£ 


^C* ^c 

0, p, r 


r 


R 


H 3 


e 

r 


angle between plane of rotation and center line of blade 
(palled flapping angle), positive upward, that is, 
from x to z 

angle between x— axis and projection of center line of 
blade on plane of rotation (called lagging angle), 
positive backward, that is, from y to x 

values of 0, 0, £ in steady state of flight 

values of 0, 0, £ if steady state of flight is hovering 

state (subscript c = Constant) 

deviations of 0, 0, £ from values in steady state of 

equilibrium, that is, from 0 O , 0 O , £ o 

distance of blade— hinge center from axis of rotation 
(see fig. 1) 

distance of point of center line of blade from hinge center 
(see fig. 1) 

tip radius of blade 


_ - ° 

% e 


v x ,v y ,v 


R e 

r i 

v 

^ = — J 


vector of resultant velocity v + TS 7 in steady flight 

, rectangular components of V 

real part of complex oscillation frequency p 

inner radius of blade length (see fig. 1) 

also a geometric constraint constant in case C, OSC TT.T.A T'TOWfi OF 
BLADE SYSTEM IN HOVERING, and IN LOW-SPEED TRAVELING 




e 



that is, speed ratio 



Si 


imaginary part of complex oscillation frequency p 
defined by equation (13), 


k 


k ~ 1 
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c 

°l/ c l 

^nc-, 


°1 


- ^°S ^ac 


M xl ,M yl ,M zl 

B 

A 

V T 2 
T 12’ \ 


chord of cross section of "blade 

chord of innermost cross section (r = r A ) 

■ uni t vector of zero-lift line of cross section of blade 

component of relative inflow velocity in plane of cross 
section of "blade 


component of V, 


nc-, 


in direction of c. 


component of "V" nc ^ normal to c"^ 

circulation around cross section of "blade 

distance ratio between center of gravity and aerodynamic 
center of cross section of blade (see fig. h) 

moment components about Xj— , y-j— , z^-axes, respectively 
(see fig. 3) 

moment of momentum vector of blade 

area of blade cross section 

moments of inertia of cross section of blade 


volume integrals of moments of inertia (see equation (6l)) 


T H =l 


Ar 2 dx 


r i 

iB 


S = 


Ar dr 


P 

a 


^•complex 

^complex 
fl 


Pn 


4 n 


density of air 

density of blade material (average density in case of 
framed structure) 

complex circular frequency of oscillation 
complex frequency ratio + 1 S i,n) 

! ft \ 

real frequency, cycles per second ( Sj n , 

\ / 
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log 


*n 

Si+l 


e.n 


logarithmic decrement [ 2 it — - 


B, F, D 

Kj A, p 

E = ir § fl% 4 c i 


S i,n / 

complex amplitudes of oscillation of 6 , P, £, 
respectively, about steady state of flight 

constants in kinematic constraint conditions 

Lagrange multiplier 


* £ 5_£i 1 1 + ia 4 
0 Lg 315 V 2 1 2 0 


c 2 = 


315 


1 1 \ £ B ^°i 


v 


1 + 2tj 


k l ,k 2 ,k 3 

g 


values of integrals given by equation (21) 
acceleration due to gravity 


AERODYNAMIC AND CENTRIFTJGAL FORCES AND MOMENTS 
IN STEADY HORIZONTAL FLIGHT 


The aerodynamic forces acting on a blade when it is rotating in a 
conical path and &1bo moving horizontally will be determined basically 
by means of the Kutta-Joukowski lift theorem. This principle requires 
that the components of the total relative inflow velocity V and also 
of the circulation F first be obtained, 

A right-hand Cartesian coordinate system will be used, in which the 
z-axis coincides with the axis of rotation, or axis of the cone, and is 
directed upward, though not necessarily exactly vertical. The axis of 
the conical path of a blade is supposed to coincide with the z-axis, 
which, from an accurate standpoint, implies that the connection of the 
hub with the driving shaft is such that the blade system tilts with the 
hub. The x-exis, moreover, is in the direction of the arm (of length e) 
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of the huh. (See fig. 1.) Instead of using the direction of f light , 
which is assumed as horizontal in this paper, as the axis of reference , 1 
it is equivalent and more convenient for the kinematic analysis to con- 
sider the longitudinal axis of a blade as fixed and the line of action 
of the velocity vector as rotating. The. velocity vector will, in general, 
not lie in the plane of rotation . xy but will make an angle 7 with it. 
(See fig. 2.) Thus, while the direction of travel, that is, of v, is 
assumed to be horizontal (or nearly horizontal), the plane of rotation xy 
will be tilted at an angle 7 to the direction of flight (travel), so 
that the normal z of this plane may be tilted forward in the direction 
of flight and possibly also sideways. It is possible in this way to 
treat the motion of all blades by analyzing only one as a sample. 

As seen from figure 1, a concentric arrangement of the flapping and 
lagging hinges at distance e from the axis of rotation is assumed. Such 
an arrangement appears to the authors to have the advantage of transmitting 
the large centrifugal force on one (spherical) bearing surface with lower 
(unit) pressure than in two smaller surfaces of the sleeve bearings. More- 
over, it is shown in this analysis that the actual value of the small 
distance between the hinges has little influence on the stability 
characteristics of the rotor system. 


\ 


* 


■w 


Relative Yelocity Components 


Let v x , Vy, v z be the components of the relative (travel) velocity 

of flight, that is, relative to the x, y, z coordinate system, which is 
for this purely geometric discussion considered fixed in Bpace. Moreover, 
let v , v , v be the components of the relative rotational velocity 

of a blade. Then the components of . the resultant inflow velocity (excluding 
the induced velocity) will be 


V = Y + v 
x X fix 


T y - 


T y + T 0y 


(1) 


Y 

z 






As in or din ary wing theory, the induced flow across a helicopter blade 
will be taken into account by the changed (induced) direction of the 
total inflow; whereas the induced change of magnitude of the velocity 
is negligible. ; 
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The following notation is used in figures 1 and 2: 
ft • angular velocity of hub due to torque of- engine 
r distance along blade measured outward from hinge point H 
r centroid axis of blade 

3 flapping angle, that is, angle between line r and xy— plane 

£ lagging angle, that is, angle between x-axis and projection of r 
on xy-plane (It may be noted that the positive direction of £ 
is opposite to that of ft . ) 3 

+ angle periodically traversed by velocity vector, that is, angle ___ 
between x-axis and projection of relative velocity of flight v 
on xy-plane (in accordance with the kinematic inversion) 

7 angle between relative velocity of flight v and plane of rotation, 
that is, xy-plane 

From figure 2, the components of v are: 


v = — v cos 7 cos \jr 


v = — v cos 7 sin i|r > (2) 

v 

v z = — v sin 7 


From figure 1 and the fact that the rotational velocity Vq is perpendicular 
to~the line r^ and to the z— axis, the components of v^ are: 


v = — ftr cos 3 sin £ 


v = — fi(© + r cos 3 cos 5) ^ 


( 3 ) 


V ftz " 0 


In practical cases 3, £, and 7 will be small quantities, in the 

sense that powers of- these quantities above the second can be neglected. 
Moreover, when it is considered that the appreciable contribution to lift 
nnd torque will be made only by the blade farther out from the hinge, it 
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is permissible to treat e/r also as a small quantity for all values 
of r contributing to the forces. That is. 


P « 1, £ « 1, y « 1, tj << 1 


( 4 ) 


where 


n = 


e 

r 


Thus, by adding equations (2) and (3), making use of equation (4), and 
introducing the dimensionless variable 


fir 


the components of the resultant relative velocity, to second powers of 
small quantities, are seen to be: 




V y = - to 


i + n f 1 - 1- j c °s ♦ 




l + u 


^1 - Bin + 4. n - i + B 2 ) 


> 


V = — &r \iy 
z 


(5) 


Lift Components 

According to the Kut ta-J oukowski theorem the lift per unit length of 
a blade is given by the vector product 


L« = p T x r (6) 


where p is the density of the fluid medium, V i3 the vector of the 
resultant relative velocity (see equation (5)), and p is the circula- 
tion vector. The direction of p 1 coincides with that of the bound vortex 
line representing the blade; that is, the direction of J r is the s ame a3 
that of the vector r. The magnitude of the circulation p is given by 
the condition of finite velocity at the trailing edge of the zero— li '"t 
chord c-> of a blade. The radial component, that is, the component of 
the velocity parallel to the blade axis r, will have no influence 
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on the value of F, "because the circulation Integral following this 
radial component, which is the same above and below the blade, is zero. 
Thus F itself is determined by the velocity component in the plane 
perpendicular to the blade axis; whereas the lift (see equation ( 6 )) is 
determined by the total velocity vector (including the radial component). 
A well-known formula gives the magnitude of F by 


r 


jreV 


nc 


a 


( 7 ) 


where 

c chord of cross section 



velocity component of V in plane perpendicular to r, that is, 
in plane of cross section of blade 


a angle of attack, that is, angle between V' and c^, line of 
of zero lift of cross section 1 


The velocity component V 


nc-, 


can be determined as resultant of 


V n and 


and V, 


C 1 J 


where 


■n 


component of V in plane of blade cross section (at any r) perpen- 
dicular to zero— lift line 


V c ^ component of V in plane of cross section parallel to c^ 


If 0 denotes the pitch angle (that is, the angle between the zero— lift 
line of an airfoil section and the plane of rotation), which like P, t, 

7, and 11 may be considered a first— order small quantity, then, with 
the direction cosines of the vectors and r (see appendix A), the 

expressions for V c -^ and V n , to second— order small quantities, are found 
to be: 


V Cl = fir 


1 +■ u (l — - — I sin t + T| 


2 


+ u(£ — 00) cos ^ + \xQy — 


P 2 + 6 2 


’n 


Hr-j^tl + t|) +■ u[(p + £ 9 ) cos V + (0 — 3 £) sin ^ 


(8a) 

(8b) 


* 


9 


r 


9 


The angle a can be obtained from the relation 
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a 


tan a 



( 8 c) 


whereas 



can he obtained from 



From equations ( 8 a) and ( 8 b) it can be seen that V n is first order 
small , whereas V Ci is finite j it follows from equation ( 7 ) therefore 

that to secondr-order small quantities 


r * rtcv„ 


(9a) 


or, from equation ( 8 b), 


r = itcftr <0(1 + ti) + d [jp + £0) cos ♦ + ( 0 — P£) sin ♦ — r} ' (9b) 

Therefore, by expanding the vector product of equation ( 6 ), using 
the velocity components (equation ( 5 )) and the direction cosines of r, 
the lift components per unit length of a blade, to second orders, are 
found to be (see appendix B): 


V 

y 


V 


— itpcft 2 r 2 p(l + H sin +)[0 + d(P 00s + + 0 sin \|r — 7)} 

— 3tpc^r 2 u(7 — P cos if) [0 + d(cos + + 0 sin ^ — 7 )j 

nc^r 2 p'< 0(1 + 2r j) + d[(P + 2 £©) cos ^ + (20 — p£) sin t — 7J ^ (10 

_ 

+ t|d [p cos ^ + 20 sin if — -0 + d 2 (sin if + £ cos if) 

x [©(sin ^ + 5 cos ♦) + p(cos if — £ sin t) — 

It will be observed that, as might be expected, 1 ^* and L • are 
of a higher (second) order small than L z * (which is first order small). 
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Drag Components 

Tlte drag components per unit length of span of a "blade can "be de- 
termined as follows: Inasmuch as the drag D* will be parallel to the 

resultant relative velocity V, it follows that 

D ' « D*I^ D • . D*2e, D * » D*~ (11) 

r V 7 y z y 

The total drag per unit length can be expressed by the equation 

D * " ^'(“p + a i) (12) 

where ctp and are the parasite and induced changes of the angle of 

attack, respectively. It may be remarked that equation (12), although 
not exact, is correct to quantities of second order, since, by observing 
the second-order smallness of L^.* and 1^* , the total lift to quanti- 
ties of second order may be given only by the z-component L z *. 

From equations (11), (12), and (5) it is seen, by considering 

and, oj as first order small, that D z * will be a third-order small 
quantity. Hence to second orders. 


V - 0 


As can be seen from equation (12), it is sufficient, in order to determine 
D x * and Dy* to second orders, to consider only the finite terms of 

V x , Vy, and V" z . Thus, from equation (5), 


* \ Jv x 2 + V 2 = Hr(l + d sin ijr). 1 + ( 008 ^ ^ 
' \l \l + u sin^y 


For purposes of investigating stability it will suffice, in order to 
avoid needlessly complicated integrations, to replace the radical factor 
in the foregoing expression for V by an average value. Thus, taking 


u 



as the highest expected speed ratio. 


cos ~ sin^V ~ i 
2 * 2 
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it is seen that the Itiaximum. of 


p. cos jr N 
1 + y. sin ty 


1 1 + 0.11 ~ 1.087 


for the highest expected velocity vf^ = 2 / 3 ). Therefore, there can he 
written ' 

V ~ &r(l + \i' sin ijr)fc (13) 

where k varies from 1 (for hovering, or v = 0 ) to I.O 87 ^for high-speed 

ratio — = nj) . Thus, hy use of equations (lQa), (ll), (12)* and (13)* 
flR e / 

the expressions for the x— and y-components of the drag per unit length, 
to second-order small terms, are found to he: 

D x * m — — (otp + Ojjc^r^pp. cos + ju(0 cos ^ + 0 sin ijr — 7) + ©j 

i. U l 

D y * = - ^(ctp + a i ^cfl 2 r 2 p(l + n sin *)|u (0 cos ljr + 0 sin + - 7) + 0 ] 

For purposes of investigating stability, when the assumptions or 
approximations need not he as accurate as for performance, the induced 
angle of attack may he approximated hy the average along the radius, 

which can he taken from the well-known formula for elliptical wings of 
aspect ratio AE in rectilinear flight: 

2jc 




U L = AE = 2 

rt AE « AR AR + 2 


The angle of attack a can he determined from 


Hence, hy using equations ( 8 a) and ( 8 b) and supposing AE = 8 , the 
expression for the induced angle, to first orders, becomes 

i 0 + jj.(g cos ^ + 9 sin ijr — 7 ) 
a, = i 

' 1 + u sin ^ 


( 15 ) 
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When, for simplicity, the case of hovering (m.« = 0) is considered., 
it will he noted from equation (15) that if the blade angle 9 is assumed 
as constant along the blade (see next section), then the induced angle of 
attack has implicitly been assumed to be constant. However, because the 
chord will be taken as a parabolic function along the blade decreasing 
toward the tip (see equation (l 8 )), the drag force per unit length will 
also decrease parabolically. According to the theory of minimum drag of 
rotating airfoils, the induced angle will increase slightly near* the root 
and then decrease, also slightly, toward the tip, according to a function 


of the form proportional to 



where w is the axial com- 


ponent of infxow velocity, which is equal to the induced velocity in 
hovering. When the multiplication by the chord is considered, a more 
exact determination of the drag force according to the aforementioned 
induced inflow distribution cannot cause much difference from the results 
of the assumptions made in this paper, particularly because for the blade 
treated as a rigid body oscillating about its hinge, only the resultant 
moment of drag enters the dynamic equations. 


Pitch Angle 

The blade angle is determined by the weight which the blade has to 
carry by means of its lift force. For this purpose it will be sufficient 
to consider only quantities up to first order small. Thus, from equa- 
tion ( 10 a) , 


V * 0 

L * = 0 

y 


V 


r — 

itpcfi 2 r 2 (l + u* sin 10 6 + n(f3 cos if + 9 sin If — 7 ) -fi^aAr 


d g 3 f 
3* 2 


( 10 b) 


The last term of equation (10b), where A is the area of a cross section, 
represents the inertia force for the case in which, for p e / 0 , the 

flapping angle 3 0 in steady flight varies periodically with if( = nt), 
that - is, with time. The ideal requirement would be that- the lift of a 
blade element be constant during a cycle. The pitch angle. 9 would then, 
according to equation ( 10 b), have to vary (with 10 along the circumference 
as follows: 



NACA TJT No. lU-30 


15 


0 


hi + a S 2 Ar 



ii(P oos t - r) 


(Hr ) 2 * 0(1 + » sin *) 2 1 + I* =ln * 


( 16 ) 


In equation (16) L z *, £lr, p., P, and y must be considered as 

constant on the circle at radius r vhile ijr is changing from 0 to 2jt. 
■Equation (l 6 ) would require a freely twisting blade, which, of course, 
would be difficult to realize. Therefore, it may henceforth be required 
that only the total force component tf L z * dr be constant (t hat is, be 

independent of *). Then 0, which will be dependent on *, can be 
determined as follows: Assume 


0 »» 0 C + AS 


(17) 


where 0 C = f(r) and A 0 = g(i|r, p, 7 , v) bub is independent of r. 
This means that the blade is designed sufficiently rigid so that, for 
change of pitch, it can practically be rigidly rotated only about its 
radial axis. 

Let rj_ be the radius of the Innermost section of the blade, or 

the value of r at which the blade begins to beccme effective in lift. 
The following dimensionless quantities, moreover, in which R is the 
tip radius of a blade, will henceforth be used: 



As an example , it will be assumed that the variation of the blade chord c 
with the radius r is parabolic. This provides 'a decrease of the lift 
to zero at the tip of the blade and represents in a way the decrease of 
circulation toward the tip, without making it necessary to enter into the 
theory of the trailing distorted helical vortex sheet. Thus, it is supposed 
that 
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c = 



(18) 


It seems appropriate, moreover, to choose for an example the "basic pitch 
angle 0 C for the hovering state (h = 0 ) as constant (with r) . It 

should, he remarked, however, that this assumption, as well as that of 
equation (l8) , is suggested only hy way of an example to fix the ideas 
and the order of magnitude and approximate form, of the functions appearing. 
Whether these or other exemplifying reasonable assumptions are chosen will 
have only a negligible influence on the frequencies apd on the stability 
conditions. 


Put 6 = 0 C , u = 0 , and equation ( 18 ) into equation ( 10 b), and set 



L z * dr 



cos 7 


W 

n 


( 19 ) 


where W is the weight of the helicopter, end n is the number of blades. 
Then by evaluation of the integral in equation (19) , with the substitution 
of equation (,10b) , the expression obtained for 0 C , to first powers of s A , 
is found to be 


9 


e 


c 


W 109 

pR^ 8njr(2 + s i ) 


( 20 ) 


In order to find £B , pub equation (17) into equation ( 10 b) 
require the condition of equation (19) , substituting for W 
0 C by means of equation ( 20 ). Then the expression for £0 , 

of Sj higher than the first are neglected, is found to be 

” <i epo( s l n ^^ ^ 2 kg + n e k 1 siny^+( Pcoety - 7)^2 + ^e k l sini M 
£0 = : Z 


and again 
in terms of 
if powers 

+ £ _S_ 

P ^ 


1 + 2u e k2 sin* ' + fi Q 2 k-|_ sin 2 i|f 


( 21 ) 


where 
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With the assumption (18), the values of k-^, kg, and k^, to first 
povers of Sj_, are: 



Only the total lift force L z has "been calculated for the purpose of 
determining the "blade angle. The other force resultants L , L , D . 
and Dy follow, "by integration, from their unit values (equations (10a) 
and (14)), and determine, together with L z , the reaction of the hub 

and the fuselage. They are not needed in this paper, which confines 
itself to the blade system. 


Aerodynamic and Centrifugal Moments 


In order to determine the steady— state flapping and lagging angles 
P Q and and to investigate the stability of the helicopter rotor, 

it is necessary to obtain the moments acting on the blades. The moments 
about the following axes must be determined (see fig. 3): 

The x^-axis, giving the twisting moment which in some way must be 
held in balance by the pitch— changing mechanism of the hinge 

The yj— axis (perpendicular to the Xj— axis, and, like x^, parallel 
to the xy— plane), giving a condition for the flapping angle p Q 

The z-axis, determining the torque per blade 

The z^-exis parallel to z and passing through the hinge point, 
giving a condition for the lagging angle £ Q . 

Twisting moment about x^-axis.— By referring to figures 3 and 4, it 
is found that: 


M xl - **ac - 






L* cos 9 cos p 
« 


dr — afP 



I 


P 


- (* C g “ *ac) ( L xl* + D xl') 008 9 sin 9 


dr (22a) 
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where M ac Is the (aerodynamic) moment about the aerodynamic center and 
can he expressed as 

M ac -/r" ^ t*! 2 ' 2 

The density of the blade material is 0; and I ? are the principal 

moments of inertia of a blade cross section ^1^ > I As is well known, 

i (i-j^ — 1^ sin 2 0 is the product of inertia about the x g — and y 0 -exee 

(referred to in equation (60)) of a cross section. In the following 
text the abbreviation *12 - J 1 dr ^I 1 — will be used; Ip denotes 

the polar moment of inertia of a cross section about the y B — and z 8 -axes 
(see equation (6l)); I cg and l (see fig. 3) are the distances of 

the centroid and of the aerodynamic center, respectively, of a blade 
cross section from the leading edge of the airfoil section. The abbrevi— 
l - 1 

ation f s ■■ - a - C will henceforth be used. The last term of equa— 

c 

tion (22) represents the twisting moment of the inertia force for the 
case in which, for u© / 0, the pitch angle 0, in steady flight, 

varies periodically with Ht), that is, with time. Analogous 

inertia terms will also appear in the expression for the moments M yl 

and M zl , equations (23) and (29), respectively. These inertia terms 

are all derived in detail voider INERTIA FORCES AND MOMENTS AND EQUATIONS 
OF 0SC U. I ATION and are given in advance here only to make the expressions 
for the moments M^, M^, and M z i complete. 

The second term of equation (22a), represents the moment contribution 
of the distribution of the centrifugal forces in the cross sections. The 
third term of equation (22a) is due to the fact that for rigid— body oscilla- 
tions the centroid axis (see fig. 3) (which will generally not pass through 
the aerodynamic center) must be taken as the reference axis. For elastic 
vibrations with the twist angle equal to zero at the root the shear center 
(which for airfoil sections, however, lies very near the centroid) must 
be substituted for the centroid. 


By keeping in mind the orders of magnitude of the terms involved and 
rejecting all terms smaller than the second order, remembering that 
Cw will, in general, be negative, and assuming (as appears reasonable) 

( l cet “ l ac\ 

that \ — = f is constant ( ~ 0.15) along the blade length r. 


the expression (equation (22a)) for may be written as follows: 
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T Cl 2 c 2 dr + f I L z *c dr 


-ati 2 e 


f* (i, . ^ ^ 0 /; 


I p dr (22b) 


Moment about y-j-axis It will be seen from figure 3 that this moment 


-r tv r cos P — * + D ri *)r sin pj dr 


+| r P cos 


(5 - 5^ ao * on 2 |4 x s 


where % H / Ar dr. 


In equation (23), the differential dC of the centrifugal force Is 
given by 

dC = n 2 ^ dm = n 2 r L A cr dr (2h) 

where A is the cross-sectional area of a blade element at the distance r. 
The value of r^ can be obtained from 

r^ 2 « e 2 + (r cos P) 2 + 2e(r cos P) cos £ 


% r(l + n - | 3 2 j 


to second-order small quantities. 
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The value of ~ Zj) C£m b© determined (see fig. 5 ) from the 
approximate relation 

(* 

or 


- Qrs; (?r) 


e 


e + r 


(f “ *1) ~ 


( 26 ) 


to second— order small quantities. 


When equations (24) , (25) , and. ( 26 ) are put into equation (23a) and 
terms smaller than the second order are rejected, the expression for M 
"becomes y 



/ _\ d 2 3 _ _ 

dr + + cQ 2 — I H + agS 


(27) 


where 


S Ar dr 

U r i 


The last term of equation ( 27 ) represents the effect of the weight of a "blade. 
Moment about z^— axis.— From figure 3 , this moment is seen to "be given "by 

M zl (^l* + D ylj) r to* +J r sln(s - &^dC + Ig ( 28 a) 


From equations (24) and (2 6 ), to second-order small terms can "be 

written as 

M zl “J r ( L yl* + D yl') r ** + S + an 2 ~ 




( 28 b) 
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Finally, the moiasnt about the z— axis will he 
nE 




0 R 

(jyi* + D yl‘) r l 003 (f - 0 dr (l^ + r x sin (S - ^ dr 

(29a) 


S 2 t 


To second-order small terms, this may he written simply as 

>R 


M„ 


i *>' + V) 


r dr * M zl — oJi 2 ^e S 


( 29 b) 


It may he remarked that the center of pressure on a blade at any speed 
ratio M. e moves radially with the angle of rotation V. This motion of 
the center of pressure is, however, taken into account by the moment- 
equilibrium. relations between the aerodynamic, centrifugal, weight, and 
inertia forces. This radial displacement is compensated for by the change 
of pitch angle Q with the azimuth angle which can be accomplished 

by either a swash plate (sufficient near hovering for the first harmonic 
term) or a cam plate (for higher speed ratios). At higher speed ratios 
(not treated in this report) the problem of the blackout of lift toward 
the root of the retreating blade would have to be considered. 

Figure 3 shows that any force components in the x^— and y-j— directions 
are related to those in the x— and y— directions by 


F xl* * V 008 £ - F y* sin 5 , F yl * - F ' cos £ + F x * sin C 


When it is remembered that F * and F * 

i J 

small quantities L^-, 1 , ^yi* •> , and 

small, it follows that to second orders 


are in this case the second-order 
D* and that £ is first order 

v 


F 


xl 


t ~ 


« F. 


F yl* ~ 


F * 

7 


This simplification will be used in the following section. 

Explicit expressions for moments .— The four moments, as given by 
equations (22b), (23b), (28b), and (29b) can be evaluated by using 
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equations ( 10a) and (ik) for the lift and drag components and the assump- 

tion (i8). for the values of the zero-lift chord c in terms of rfor s) 

that is, along the "blade axis. In the evaluation of these moments, for * 
the analysis given here, all terms smaller than the second order have 
been rejected, and for further simplification, even second-order terms 
have been neglected when they appear as additions to first-order or finite 
terms (for- example, as in L z * ) . Moreover, p overs of sj higher than 

the first have for simplic±ty been neglected. The results obtained for 
the moments are then (see appendix C): 


M xl " -P fl2R3 °n 


Hnfl E- + s i + ^e (l 


si) sin y + 3 ^ e 2 ( 1 - s i)(l — COB 2+ )j 


•■.>[&• .4 


9 sin if - 


0 

+ -jj— (1 - 2s i )(0 sin 2y + 0(1 - cob 2i|r) - 2 y sin if) 


y^9l - CT n 2 Ar, dr 

12 ^2 d P 


l + f 008 *) 

3 


& 2 u P 


( 30 a) 


For hovering (n 0 a 0 , 0*0), 


W. ’ - (1 * [|°«J - 


-ail 2 0_ I 


c - 12 


( 30 b) 


where, as before: 


*X2 m 


r & - *«) 


dr 


M. 


yi 


315 


^1 + ^s^rtpO 2 R^c^|80 + 12n e (p cos if + 20 sin y — 


+ 21u 0 "(p oos y + 0 sin y — 7) sin yj 
+ + Sej + an 2 — ^ X H + Sag 


( 31 a) 


ay* 
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For hovering. 


^My.^ » - itpSpR^c^ 0 O + oft 2 P + Sog (3113) 


M zl “ “ i 


8e 

315 k 


(«* + I 9) + i*. ^ [ s (r - * co« * + £ Bin *) 

+ 5^5 + y s )(0 cob ♦ + e sin * - rij - H 1 * 2 008 * 

— j 

+ 0 sin * - 7) {7 - P cos i ~ Bin * + 008 * + 9 sin ♦ " r ) 


+ ae + otfi Xg 
Si * 2 


( 32 a) 


For hovering. 


( M zi) c * - 35 ” 0 *S «c(^ * \ e c) ♦ 5 <»> 


M z 50 Mzl “ <*&£e s 


(33) 


Expressions (31a) and (32a) for the moments M yl and M^, respectively 

will he used in the section STEADY STATE IN HOVERING- AND IN LCW-SEEED 
TRAVELING to determine the steady-state values of (3 and £, which will 
he constant in hovering, and functions of if and u e in traveling. These 

expressions for the moments will also he used in the following section to 
determine the "quasi-elastic" moments (that is, moments depending on the 
deviation from the steady-state equilibrium position) during an oscilla- 
tion. 


FORCES AND MOMENTS DUE TO SMALL OSCILLATORY 
DISPLACEMENTS AND TO VELOCITIES 


In order to treat fully the questions of f: *i ,S?i£££nt s 
and stability, the forces and moments due to oscilla y P 
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A3 s 3, AC s f , A© 5 ? and velocities P © a 

do at dt 

mist be determined. These will be the quasi-elastic and the aerodynamic— 
damping terms in the final dynamic equations of oscillations about a 
state of steady motion. The oscillatory displacements and velocities _ 
which must be considered herein are the angular displacements 3, 9 

and velocities 3, £, ©. Because only natural (or free) oscillations 

about a state of steady flight are treated in this analysis, the fluctua- 
tions (in magnitude or direction) of the rotational speed as caused by 
change of engine torque and of the translational inflow velocity v Q 

as caused by gusts or change of angular position of the driving shaft 
will be left for a later research. These latter fluctuations would give 
rise to forced oscillations, which could be treated by the same general 
method as given herein for the free oscillations. 

In the following analysis, squares of oscillatory displacements and 
velocities will in all cases be neglected. 


Quasi-Elastic Moments 

The quasi-elastic terms, which must be used in the dynamic equations 
of oscillation, can be determined from expressions ( 30 a), ( 31 a), and ( 32 a) 
for the moments by putting 


0 » 0 O + 0, 3 = 3 0 + F, £ “ £ 0 + T (3M 


in theee expressions. In equation 3^* the subscript o denotes the value 
for the steady state, which is constant for hovering (u e = 0) but a 

function of ♦ for traveling (t e 0), whereas the bar denotes the 

small oscillatory changes (varying with time) about the steady state. 

The quasi-elastic terms are the moments due to the changes 0, 3, and 

By putting, therefore, equation (3^) into equations (30a) to (32a), 
subtracting in each case the moment for the steady state (0 = 6 Q } B ■ 3 0 ; 
and £ ® £ 0 )> and rieglecting powers of the oscillatory changes higher 
than the first, the following expressions are obtained for the (aero— _ 
dynamic and centrifugal) quasi-elastic terras: 


0 + 2u e (20 sin t + 3 cos t) 


(V) 




qe 


12 


f 1 


+ 


3u e 2 (l — 2sj) & sin 2!f + 0(1 — cos 



V I 


12 


(35) 
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cos if + 20 sin 


*) 


+ 21 cos + s’ sin if) sin ty] + er^ 2 ]? ^1^ + S ej 


(36) 


( M Z l\ e " " ~ 0 (a p + § 0 O ) + ^e 155 [° ( 7 " 


B 0 cos if 


+ ^ sin i^ - F e o + |(F cos if + 0 sin y)^a p + | e^j + 0 (3 

+ 0 O sin if - r)J + t^ e 2 ^ |(F cos if + 0 sin +) j^(p o cos if 


cos 


CL 

+ 9 Q sin if — r) + y - P Q cos t + ~ sin i|r 


+ 0 Q sin t - r) [■ ] + cr©n 2 ^ S 


— (3 cos t)(P 0 008 ^ 


(37) 


Components of Total Be lat ive Inflow Yeloc ity of Oscillation 


By referring to figure 1 and using the fact that, at any distance r 
along a "blade from the hinge, the linear velocity rP due tojangular 
oscillations of the flapping angle will he perpendicular to r and will 
lie in the plane which contains r and is perpendicular to the xy-plane, 
it is seen that the velocity components due to rP will he: 


Y x £ » r 3 sin 3 cos £ Z r 33 




z3 


r 3 sin 3 sin £ » 0 r 


r 3 cos 3 ft; r 3 


(38) 


S imi larly, use of the fact that the linear velocity r£, due to oscilla- 
tions of the lagging angle, will he parallel to the xy-plane and perpen- 
dicular to the projection (of length r cos 3) of the hlade axis r on 
the xy—plane yields for the velocity components of r£: 
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r cos 3 sin £ * r^^l 



r cos 3 cos 


£ * *£ 


(39) 



J 


All approximations mads in this entire section are up to second— order 
small terms. 


Damping Forces 

The contributions to the aerodynamic forces following from the 
changes of the velocities calculated in equations ( 38 ) and ( 39 ) can he 
determined in the following general manner. The velocity increments 
must he added to the corresponding expressions (equations ( 5 )) for the 
velocity components. The forces and moments are then determined hy the 
same procedure as in the section AERODYNAMIC AND CENTRIFIXiAL FORCES AND 
MOMENTS IN STEADY HORIZONTAL FLIGHT < and the damping terms will he those 
containing the quantities £ and 3 in the expressions for the forces 
and moments. In this manner (see appendix D), the expression for the 
total circulation to second orders is found to he: 


T s Kcflrje(l + t|)+u(3+S0) cos ♦ + (9 — 3£) sin \|r — - 0 




Q Q 


(40) 


Then, hy the use of the Kutta-Joufcowski relation given in equation (6) 

(see appendix D), the increments due to 3 and £, in the lift components 
per unit length, are found to he 




AL* ■ 
z 


AL = itpG 2 r 2 o3(l + U sin >10 — 

. n 

0 

— . rtofi 2 r 2 c ■§ + m( 23 cos \|r + 9 sin ¥ — 2 y)J 

% 

— ffoG^^c 4 (1 + 71 + p£ cos 7 ^ + [i. sin ♦) &■ 

n 

+ [2 9 + u(3 cos t + 9 sin ♦ — 7 )] 1 

0 " - 1 


(4la) 

(4lh) 


(blc) 


The additional drag components per unit of hlade span, to second orders, 
can in accordance with equation (ll) he determined from 
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AD * 


AT\ t- 




(42) 


where £D* is the increment in the resultant drag and. may be given, 
according to equation (12), by 


AD* - (c^ + ct ± ) AL Z * 


(43) 


Therefore, by putting equation (4lc) into equation (43) and then putting 
equations (43), ( 5 ), and ( 13 ) into equations (42), the expressions for 
the additional drag components are found to be: 


AD. 




jtcn 2 r 2 M.(cos +) ^ 


3 


AD 


AD Z ' 


' • - i l «c pn^d + 

y v 


U sin +) ' 


(44) 


As in the section AERODYNAMIC AND CENTRIFUGAL FORCES AND MOMENTS IN STEADY 
HORIZONTAL FLIGHT, may be obtained from equation ( 15 ) . 


Damping Moments _ 

Damping moments about the hinge point are caused by damping forces 
distributed over the length of the span of a blade and acting with their 
radial levers. Such moments appear furthermore as the effect of the 
pitch-angle oscillation. This latter effect, for the moderate velocities 
appearing in helicopter flight, can be calculated by means of an apparent 
change of local angle of attack under the assumption of quasi— stationary 
flow. This calculation (see appendix E and fig. 6) leads to the relation 



where, from equations (8a) and (l8), 
>R 
1 r 


r dr 7 ci c 3 = r ® 2 ^ 3 ~ (l + s^ [2 + s 1 + ?n e + 7 U e (sin * + Z cos V)j ( 45 ) 

J r l . 
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This simplified approach is made for the damping moments, although it is 
well known that for high values of osoillatory and inflow velocities an 
effective change of camber of the airfoil, as well as the reaction due 
to the trailing vortices of the flow around the airfoil, would have to 
"be considered. As has "been made plausible from the theory of straight- 
moving airfoils, however, these effects are small for velocity changes 
occurring in helicopters. A sufficiently exact theory, moreover, 
including radial velocities and helical distorted vortex sheets has 
not yet been developed. 


Another problem which must be considered in a further development 
of the dynamic theory of flexible blades, with or without fixed roots, 
is the influence of phase differences between flapping and twisting due 
to the distance between the center- of- gravity axis and the elastic axis. 
This problem does not appear with hinged and sufficiently rigid blades, 
which are assumed in this paper. 


The damping effect of flapping and lagging oan be obtained from the 
expressions for the damping forces per unit of blade length in equa- 
tions (41a), (4lb), (4lc), and (44). If the moments of these unit forces 
are integrated over the length of a blade analogously to the method of 
equations (22b), (23b), and (28b) for steady flight, the following 
damping moments, including also the effect of pitohr-angle oscillation, 
are obtained. The increment in due to damping will be 








dr + f 


* T AL * c dr - — p0 r V 1 

J-i Z 12 Jr i Cl 


c? dr 
(46) 


where Y, 


C 1 


AY, 


AY„ 


'x* ^"y* 

cosines l 


(see appendix F)_ is to be determined from the components 
and AY Z of -equations (38) and ( 39 ) and from the direction 
and- n„_ (see appendix A). Similarly, from 


ci* m C]_* “ci 

equation- (23), the increment in is seen to be 


*Vi = -f R *v 

Jri 


r dr 


(^7) 


Finally, in accordance with equation (28b), the aerodynamic and centrifu- 
gal < 3 ampins: moment about the z^— axis will be 


am, 


zl 


= fr C^ 1 ’ + ^ Dyl *) r ** 


( 48 a) 
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As shown "before. 


V 


(or v)*V( or V) 


to second— order small quantities. Therefore equation (44) may, to second 
orders, "be written thus: 


AM 


zl 


Sr, K’ + 


r dr 


(48b) 


Putting equations (4la), (4l"b), (4lc) and (44) into equations (46), 
(47), and (48b) , neglecting terms smaller than the second order (taking 

and all oscillatory velocities as first order small), and rejecting 

powers of s i above the first results in the following explicit expressions 
for the aerodynamic damping moments (see appendix F): 


AM 


xl' 


p n% 3 Cl 2 

_ 


C 1 + s 3)( 1 71 r % & + 8 i + K 1 “ 2s i) ^ 


+ 1 


+ 7(1 — 2s^ [i e ( sin y + t cos t)j 

f it jis + 2n e (P cos \|t + sin t — y)J — J C 

L 

+ f it|l + 2 t} q + 2n e (sin t + £ cos t)j 
AMyj_ = itpfl 2 !^^ |l + ^ [_4S + 3M- e (P cos t + 9 sin t - 7)] 


(1 + 2(i Q sin *)f 


(49) 




+ ^ £2 + 3ri e + 3u e (sin t + £ cos ♦j]f 


(50) 


AM zl = — 


+ 3H- 


— (l + - e}j I J2/1 _ 0 - 1 a 

315 \ 2 y a \ 5 V k i 

— 

- — ) 0 00s ♦ + (1 - — } 

5ty \ 5k/ 


a, / ’ 

9 sin ^ sin ty — 2 (1 — — - 

k V 5k 


K5D 
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In addition to the aerodynamic damping moments given in equations (t-9) 
to (51), there may, in general, he "pseudo-damping" moments due to a change 
in the angular velocity a caused by £. This cha n ge can be determined 
by substituting (ft — 5) for ft in expressions (30a) to (32a) for the 
aerodynamic and centrifugal moments, subtracting in each case the moment 
when £ = 0, and neglecting squares of g. It will be observed from 
equations (30a) to^(32a) that each of the moments in the steady state is 
proportional to that is 


M = Gfi 2 (52) 

where G is a function of 9 , 3, £, and so forth, but not of ft. 

Therefore the pseudo-damping moments Ap^ M will be of the form 

Ap d M = — 2G ftC (53) 

In steady flight, however, the angles £ 0 and 3 Q will so adjust them- 
selves that the moments and are both zero. Moreover, it may 

be assumed that by some means or other, for example, counterweights, the 
moment will also be balanced. Therefore the value of G as defined 

by equation (52) will be zero. It follows then from equation (53) that 
(because G is obviously the same function in equations (52) and (53)) 


^d **xl ■ ^pd »Vl " V M zl * 0 


* 

Hence, the effect of £ 


on ft need not be considered. 


(5*0 


INERTIA. FORCES AND MOMENTS AND EQUATIONS OF OSCILLATION 


The derivation of the inertia forces and moments may be based on the 
time rate of change of- the moment of momentum vector S’ referred to a 
coordinate system rotating about the driving shaft of the blade system 
with angular velocity -ft. It is well known that with such a coordinate 
system the centrifugal forces, the centrifugal moments about the blade 
axis, and the Coriolis forces must be added to the other impressed forces 
(in this case, aerodynamic and gravity forces). 

The blade system is assumed to consist of blades, each hinged to the 
hub driven by the engine shaft, and the hinge system for each blade is 
assumed to have a common hinge center, which is then the natural moment 
center for the moment of momentum vector. 
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It is true that for reasons of detail design it is often found that 
the flapping hinge and the lagging hinge are not concentric; as only the 
eccentricity of the lagging hinge is of importance, however, it is a very 
small loss of generality "but a great advantage in simplicity to assume 
concentric arrangement of the hinges. It may even, for the development 
of larger systems, he convenient to support the entire centrifugal forces 
hy one spherical hearing surface instead of hy several necessarily emaller 
cylindrical hearings, each of which must support the entire centrifugal 
force . 


Moment of Momentum. Vector 


The basic dynamic theorem, is expressed in vector notation hy the 
equation 


B * M 


(55a) 


where 


B 


is the moment of momentum vector. 


B 


aB 

dt 


and M is the moment vector of centrifugal, relative acceleration, 
Coriolis, aerodynamic, and gravity forces. 


In Cartesian coordinates, equation (55a) can he expressed hy 

K K “ V * M z (55b) 

The vector B and its rectangular components must now he developed in 
terms of the positions and velocities of the mass particles of the blades, 
whereby it is sufficient to consider one sample blade. 


The general relation in vector notation is given hy 


B 



X r 


( 56 a) 


where V-^ denotes the resultant velocity of a blade element, and r the 
radius vector from the moment center (hinge center) to the particle. 


With u, v, and w as Cartesian components of and x, y. 


and z as those of r, equation ( 56 a) is equivalent to the following 
three component equations: 


B x = / dm(wy — vz) 
B = J dm(uz — wx) 

J 

B z - f dm(vx — uy) 


(5^h) 
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The integrals in equations (56a) and (56b) must "be taken over all mass 
elements of a blade. 

The velocity vector V^, or its rectangular components, muBt now 

be expressed by the condition that the blade is kinematic ally constrained 
to move about the hinge center. This fact, in vector notation, is 
expressed by 


y-jj = r X m 


(57a) 




of the 


where 00 is the vector, with components a^., ay, and 

angular velocity of the blade. In Cartesian notation equation (57a) is 
expressed by 


u 


- y 


v = xm z — zoy 


v m yty - xoy 


(57b) 


Inserting equation (57b) into equation (56b) and factoring out the density 
or of the material of the blade gives 


B_ 


B. 


B„ 


( 5 


W 


^XZ^Z 


A - X yz-z 


- I 


(Vs - - V”y) ® 


(58) 


where the moments and products of inertia are defined (in the usual way) by 
T x => J J dA dr^y 2 + z 2 ), I y ^ J * dA dr(z 2 + x 2 ^, I z = dA dr^x 2 + y 2 )] 


yz 


II 


dA dr yz, I 


zx 


.jy’dAdr ZX, 


xy 


: 2 ), I z = / / dA dr ( x c 
xy 


I 


and where the bar over I is intended to emphasize thejvolume integrals 
as opposed to area integrals appearing later (cf. equation (6 l)). The 
coordinate system x, y, z (see fig. l) in general does not coincide 
with the centroid coordinate system fixed in the blade, except when the 
parameters 3 and £ are zero. It 1 b therefore advisable to transform 
the moments and products of inertia to the coordinate system fixed in 


the blade . 
cides with 


This transformation is simplified by the fact that x_ coin— 

s 


with c , 


') yg "J.V11 V., 

the first powers of the angles 


and 
3 


s 

and 


is normal to c and that only 
£ need be considered. Thus 
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r =» r, y « y s — £r, z = z s + 3r (60) 

The expressions (59) then take the following forms: 



and where A denotes the cross section (varying with r) of a "blade. In 
equation (6l) it has "been observed that the static moments about the 
centroid axes are zero. 


Time Rate of Change of Moment of Momentum Vector 


The tine rate of change of the moment of momentum vector B, with 
the condition that squares and products of angular accelerations and angular 
velocities shall be neglected, whereas the products of accelerations and 
the stationary angles 8 Q , 3 Q , and shall be retained, can now from 

equations (58) and (6l) be expressed as follows: 



^0%) * 

” Vo f 12 + £ 0 % 




- - a V 0 o I 12 


(62) 


The blade angles 8, 3 , and 5 111 this equation are provided with the 

subscript o in order to emphasize that they are parameters of the 
steady state, which are constant in hovering and functions of V( = &t) 
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in traveling. In order to express the moment of momentum vector by the 
angular accelerations 9 , 3, and £ instead of, as "before, by the 

angular accelerations ca z , a^, and toy, the following transformation 
may he applied by the projection of the components o^, = 9 
and = — P on the axes x and y. 


• • 

cos P cos £ + to yl Bln £ * 9 — p£ 
a^. =« — ojj, cos 3 sin £ + aj yl cos £ * — 0 £ — 3 ^ ( 63 ) 


«k--c 


The c orre sponding components of the vector B can be derived similarly 
by the relations 




Vo 


+ 



( 64 ) 


Applying equations ( 62 ) and ( 63 ) yields 



The centrifugal (inertia) forces and moments are taken into account in the 
sections giving the conditions of equilibrium of— aerodynamic and centrifu- 
gal forces, especially in steady flight (see STEADY STATE IN HOVERING 
AND IN LOW— SPEED TRAVELING) . 
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Coriolis Mcaaent Vector 

The vector a„ of the Coriolis acceleration is given "by 

a c * 2fiX V rel (66a) 

In the case considered herein the vector of rotati onal velocity coincides 
with the z— axis, that is, 

n - . 

The vector product of equation (66a) therefore consists only of the two 
components: 



hi|® 
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By use of "the transformations (equations (60)) of the moments of inertia 
to the Blade axes and the notations (equations ( 6 l) ) of the moments of 
inertia of the "blade, it is found that 


M. 


cx 


M C y " ~ 2ntTl H P5 

M cz = - 2001 ^( 3 ( 0 ^ « 2001^3 




(67b) 


where, in the last terms on the right-hand sides, the relations ( 63 ) have 
been inserted. 


The transformation to the moments about the hinge axis jj_ of the 
flapping angle 3 and about the axis r of the pitch angle 9 is again 
derived as in equation (64) and yields, to second-order small q uan tities, 
the results 


M. 


'cr 


M cy 1 “ " 2fidTg3£ 




> 


( 67 c) 


M cz = 2^33 
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The entire moment of moment tea vector is now given "by equations (65) 
and (67c), that is, by: 


B r “ 
r or 


B yl “ M cyl 


> 


B z - M cz 


( 68 ) 


Equations of Oscillation 

The equations of oscillation are formed by equating the inertia 
moments (equations (68)) to the aerodynamic, centrifugal, and damping 
moments. Thus, by using equations (65) and (67c), the oscillation equa- 
tions are seen to be: 


B r = — 3 £ 0 oI zs + fPgdjj = 

B yl - " + °^ys) + ’^0^12 + * 0 ^o^ zs + ^ P o an % * M y l + A d M yl 

B zl = - 6*(«®BC + a *zs) - 0p o a *H + P®o *12 " 2 ^o a( % - M zl + A d M zl 


( 69 ) 

( 70 ) 

( 71 ) 


In these equations £, 3 , and 5 are the total values, as given by- 

relations 0 * 0 O + O’, and so forth (equations (3^)), and the unknowns 

in the foregoing equations are the deviations 0, 3, and £ from the 

steady state of equilibrium. 

Equations (69) to ( 71 ) can be written in slightly simpler form as 
follows: From the definition of the quasi-elastic moments (given in the 

section FORCES AND MOMENTS DUE TO SMALL OSCILLATORY DISPLACEMENTS) it 
follows that, for the aerodynamic and centrifugal moments . 


M - Mqe + M Q ( 72 ) . 

where M is the total moment (that is, 9 = 9 Q + 0 , etc.), M^ Q is the 

quasi— elastic moment, and M Q is the moment in the steady state. More- 
over, for the damping moments. 
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AM « ^ M + ^ M 


(73) 


where, similarly to equation (72), AM is the total damping moment 
(9 = 0 Q + S’, etc.), ^ M le the total damping moment minus the 

damping moment Ag^... M for the steady state (0 «* 0 Q , etc.). 

Hence, using the relations (34), (72), and (73) and observing that 
the steady-state moments must by themselves add up to zero for each axis 
yields the following simplified set of equations of oscillation: 


0cl r - P^jorlzs + C0 o <tI h - ( M xl) + M . 

'dy 


(74) 


— P (dig + Clyg ] + ^0 


,) * W 0 oI 12 ♦ h 0 <ii ZB * 2Sp o o!iI H - (m,^ ♦ M yl (75) 


i («*H * 0l Z a) - + W 0 Tf 12 - SMcOSSz - (m z1 ) ?s . ^ M zl (7 6 ) 


The right-hand sides of these equations can be obtained from equations (35) 
to (37) and from equations (49) to (51). 


Relations among Geometric Constants of a Blade 

For purposes of simplicity, it will be convenient to use certain 
relations among the geometric properties (for example, moments of inertia) 
of a blade appearing in the foregoing equations and elsewhere. With the 
data of the helicopter heretofore assumed, the following relations can be 
derived ( see appendix G) : 


I ys « t m * Sr " x l 2 * s ^ 


12 x zs 


I/I H * O.OO386, so that I H + I yB » I H + I za ~ I H 

_ _ _ _ / \ R^c 

o = 2 Ij^/R, Ijj + se = Ig(l + 2*1 qJ> _ = 796(1 — s^) 

^ X H 




( 7 7j 
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STEADY STATE IB HOVERING AND IB LCW— SEEED TRAVELING 


In order to eolTe the equations of oscillation, it is necessary to 
determine first the values 0 O , 3 q , and £ Q of the pitch, flapping, 

and lagging angles, respectively, of the "blades in the steady state. 
They will "be calculated herein for the cases of hovering and low— speed 
ratio of traveling. 


If in the steady state the load-carrying condition (equation (19)) 
for 0, that is, its value 0 Q , is by Borne means enforced, then 3 

and £ , as long as no kinematic constraints are imposed on them, will 
adjust themselves freely, so that the moments and are zero. 

If kinematic constraints are imposed it is still possible, by certain 
preadjustments, to make these moments zero for the same steady— state 
values 3~ and t as without constraints. 

u o 


In traveling it has been shown that it is necessary to change the 
pitch angle 0 along the circumference of the path of a blade 
(0 ^ ^ ^ 2it) in such a way as to keep constant the force component 
perpendicular to the plane of rotation. This aim, however, will lead 
to difficulties if the blade is continued too far inward, that is, to 
r, 

a value s^ ~ so small that this inner part of the blade, when on 


the retreating side, does not find a positive relative inflow velocity. 

If this case cannot be avoided, then an additional Bideways tilting Ay 
of the axis of rotation would be necessary to counteract the moment A L_e 

Lt 

by a moment of the weight WAyh (h *» distance from, center of gravity of 
structure to center of hub). A force component transverse to the direc- 
tion of flight would arise. A loss of total lift L z would also result, 

and therefore it shall be assumed in the following discussion that L z 

can be kept constant along the circumference by means of a periodic change 
of the pitch angle 0. 


It shall not be discussed in this paper in which manner, that is, by 
which mechanism, cam device, or tilted swash plate, the variation of the 
pitch angle $ may be realized. It may be assumed that this has been 
accomplished in some way. 


Hovering 


For hovering, the steady— state values of 3 and 
obtained from the expressions ( 31 b) and (32b) for the 


£ can be readily 


moments 


and 


Ki), 




*0c 


Thus 



ko 
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Formulas ( 78 ) and (79) can be Interpreted also as an indication that it 
is possible to regulate automatically the equilibrium value of the pitch 
angle 0 C by the equilibrium values of the flapping and the lagging 
angles 0 C and £ c . 


Low-Speed Traveling 

Although the angles 0, 0, and £ will be constant in hovering, - 

they will, in traveling, have to vary periodically during each revolution 
under the action of the weight of the system and the inertia, damping, 
and Coriolis forces of the masses of the blades. It is necessary to cal- 
culate these angles because they appear as variable coefficients in the 
differential equations of oscillation. 

This calculation, however, may be simplified by assuming the speed 
ratio u e to be a quantity sufficiently small so that second and higher 
powers of it may be neglected. In the calculations it will be treated 
as a first-order small quantity. Such a simplification is advisable in 
order to gain a first access to the behavior of the blade system in the 
transition from havering to traveling. From the knowledge gained by 
such a first analysis, it will then later be possible to proceed to the 
extension to higher speed ratios. 

The steady— state pitch angle Q Q must be determined from equation (21), 
which, reduced to first powers of u e , assumes the form: 



sin ijr + 


COS T|f — 


0 



( 80 ) 



+ 



P 

OF 





where 



NACA TN No. 1430 


41 


a n d where P 0 " is supposed to he also developed up to first powers of n Q * 

(See equations (8l) and (85).) It. may he remarked that the denominator of 

equation (21) as far as it concerns the term with p 11 can he taken as 

o 

unity, since P 0 " is already itself a. second— order small quantity. (The 
damping term in P 0 * given hy equation ( 4 lc) has been neglected in this 

section hut was later found to he of the nonne gligihle order 2 . In the 
section INFLUENCE OF DAMPING ON STEADY-STATE ANGLES IN TRAVELING this 
term has heen incorporated into the 0 Q , P Q expressions. It can thus 

he finally stated that P 0 " and P Q * are the only dynamic terms which 
must appear to obtain equations (95) and (96).) 

Equation ( 80 ), however, is not sufficient to determine 0 O , because 
the <ralue of the term P 0 " appearing in it depends on the eq uil ibrium of 
moments about the y-j— axis. The moment M^ is given hy equation ( 31 a), 
which, strictly speaking, must contain in addition to the inertia term 
cx2 2 IjjP" a Coriolis term ^ C yi ■ — 2 flcIgP 0 £ 0 (see equation (67c)), a 
damping term (equation (50) ) proportional to £ o 0 q j and two inertia terms 
(equation (65)) proportional to £ Q 0 o and to 0 o £ o . Each of these addi- 
tional terms, however, is smaller than the second order and may therefore 
he neglected in equation ( 31 a). (See, however, additional term in p Q 

from equation (50) in ( 91 ).) Equation ( 31 a) for M jrl = 0 (without ( 91 )) 
gives now (hy using equation (77) ) : 


P o" + PoC 1 + ^e) 51 C 2^o + l2 ^e(fc cos ♦ + 2© c sin + - y)J - (8l) 

By eliminating P 0 ” from equations ( 80 ) and ( 8 l), the following relation 
between 0 O and P Q is obtained: 

(l + 2 n e )p o = C 2 j^20 o + 60 c + (12 - 6 k 2 )n e (p c cos t + 20 c sin t - 7 ) J ~ 

(82) 

By double differentiation of equation (82) and then substitution for P Q " 
in equation (80), a differential equation in 0 Q alone appears, name ly. 



+ 0 , 


0 C + *2^37 ~ 4 n e 0 c sin if - 2 ia e 3 c cos if 


(83) 
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The only integral of equation (83) not containing free oscillations and 
periodic in terms of sin ^ and cos V can be readily shown to be: 


0 O “ 0 c + k ^0 r “ 3 ^e 9 c sin t - | d e £ c cos if 


(84) 


This relation must therefore be assumed as the steady— state variation of 
the pitch angle 9 Q . 

If equation (84) is put into equation (82), ar. analogous expression 
is obtained for the steady-state variation of the flapping angle 3 , 
namely, 0 



Finally, the expression for must be obtained from the moment 

equilibrium about the z 1 -exis. From equation (32a) for M zl and from 

expressions (65), (67c), and (50) for the inertia, Coriolis, and damping 
moments, it is seen that up to second-order small quantities, the moment 


M 


zl 


32_ 

315 


n p ^B h c ± 0 c ^p + I 0^ 


2“ 2 *0 ~ 

+ cr£ eJTS + crn — — I H 

G*“ 


( 86 ) 


By setting M zl = 0, the particular integral (without oscillation) of 
the resulting equation is seen to be 


i 


o 




5 


c 


(87) 


Influence of Steady-State Inertia Terms 

It is interesting to observe the effect the dynamic forces and moments 
in the transition from hovering to traveling have on the values of the 
angles B Q , P Q , and £ 0 i 11 the steady state. If the dynamic terms 

were neglected, then from equations ( 80 ) and (8l) the values of 9 0 and 
0 O would be found to be (primes are used to distinguish from the correct 
values) 
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*» S- + 


c + V'e’' “ 2 Ve 8 c Bin * “ V© P c 008 * 


(88) 


IV C 1 + *>*) - Oaj^o + fas - 12 )>V + ( 12 - 


cos y + 20 c sin if 


2 g 

n% 


(89) 


When equation ( 88 ) is compared with equation (84) and it is remembered 

that with the numerical values used in this report kg = ^ it is seen 
that 0 O * is the same as 9 Q except for some slight differences in 
the coefficients of sin if and cos Thus the dynamic terms appear 
to have little effect on the value of the steady-state pitch angle. 
Comparison of equation ( 89 ) with equation ( 85 ) shows that* as for the 
pitch angle, the value of 0 O » differs from that of g Q only in the 
coefficients of sin if and cos if. These coefficients in 0 O * (without 
the effect of the inertia terms) are four— thirds times those of 0 O 
(with the inertia terms). The steady-state value of the lagging 

angle, on the other hand, is not affected by the dynamic terms (cf . 
equations ( 86 ) and ( 87 )). 


Influence of Damping on Steady-State Angles in Traveling 

In relations ( 16 ) and (21) for the pitch angle 0, the inertia 
terms of the blade appearing in the steady state of traveling have been 
taken into account; whereas the damping terms have been neglected. It 
seemed advisable, on second thought, to consider also the influence of 
these damping terms on the steady-state angles 0 Q , 0 Q , and 

(This influence, however, will not appear in the results of the oscilla- 
tion analysis for low speed ratio n e , as it would only lead to terms of 
higher than second order. For higher speed ratios, nevertheless, the 
damping terms would have an influence on the final equations of oscillation. 

Damping terms given by equation (4lc) must first be added to the 
expression (equation (10b)) for the total lift gradient L z *. In view 

of orders of magnitude, these terms may be reduced to the single term 

L * = — npfi 2 r 2 c ~ ( 90 ) 

Sty 

Addition of this term has the consequence of adding to the n ume rator of 
the expression (21) for A0 the additional term 



hk 
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(91) 


and to the expression ( 31 a) for the steady— state value of 


the term 




cM 


(92) 


Comparing expression (92) with the inertia expression appearing in equa- 
tion (31a) shows that "both are of the same order of magnitude. As has 
been seen in equation ( 8*7) ^ the value of will remain £ c to second- 

order small terms. In order to find the angle variables & Q and 3 Q 
of the steady state for n e ^ 0 (but small) the equilibrium condi- 

tions (equations (80) and (8l)) of the force L z and of the moment M 

must be complemented by the foregoing terms. Thus, instead of equations (80) 
and (8l), there is obtained 


e o - e c 




2d Q sin ^ + 3 C cos ^ — 7 J + —1— + 3 0 * 

6C 2 




3 " 
p o 


-t-(« 


~ 7 ) 




29 n sin i|r + 3 C cos f - 7 } + — - i - + 3 Q * 

8Cg 

1 + 2g 

+ + 


8c 2 


(93) 


(94-) 


Writing 


0 o - 0 c + ^e ®l» P o = P c + ^e P l 


(95) 


equating the right-hand sides of equations (93) and (9*0* and calculating, 
as before, the appropriate particular integral, it is found that the 
•value of 3 q remains the same as in equation ( 85 ), but that B Q as 

given by equation ( 8 h) must be changed by the addition of the term 
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H e Pl* = 3He C 2(3 “ 2]e 2)( 20 o cos * “ Pc sin *) 


The damping term thus appears to have an appreciable influence on 
the steady-state value 8 Q of the blade angle in forward flight, although 

it does not affect the value of the steady— state flapping angle P Q . 


OSC ILLA TIONS OF BLADE SYSTEM HT HOVERING- 


General Explicit Equations 


For the state of hovering (u = 0) the right sides of equations (74), 


(75), and (76) are appreciably simplified. 
They may be written explicitly as follows: 


2 _ •_! _ ■* _ E fl + S \ Cf Cf o / \ 

9 I r — P Solas + SPc^H + 12fi *B~ ® "if 35 \ + B * + ^ e ) 

+ jf(l + 2 0 * t(& - 5 f ^ 


(97a) 


— + Xyg j + T0 o 7l2 + ^-c^ze + 315 ^ P (l + g 6 ! + 2 le) 


+ 2T8 0 ( 1 + ! B l) + + (97b) 
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” + ^ 9 c X 12 



” 2P C X H 


+ 







(97c) 


where E a ^ £ 0 k c^. 

The external moments acting on the hinge system and transmitted through 
the hub are , Myi* and. M zl . They hare to satisfy certain conditions 

in order not to violate the equations of equilibrium. One of the possible 
sets of conditions is, for instance, W = 0(0 = 0 C ) with p and tf free. 
In this case My-, 0) gives the moment to be enforced by the pitch- 
changing lever or gear, and and are zero. 

Another condition might be suggested by such kinematic constraints 
between 0 and P, and 0 and 5 that the equilibrium positions 0 C , 

P c , and are obtained automatically without an external pitch-changing 

gear. A further condition may consist in fixing 0 by an external kine- 
matic constraint (pitch-changing mechanism) but using a kinematic constraint 
between P and t . 


It is also possible to introduce into equations (97a), (97b) , and (97c) 
friction constraints, as will be shown in one of the examples. The choice 
between these and other possibilities will be made from the following 
points of view. 


On the one hand, it is desirable to keep the natural frequencies 
away from resonance with the circular frequency SI of the drive and at 
the same time to achieve a sufficient damping decrement . On the other 
hand, it will be necessary to make sure that the kinematic conditions 
do not interfere with the transition to traveling, that is to ^ 0. 

The effect of internal kinematic conditions (constraints) will most con- 
veniently be determined by Lagrange multipliers. 


The kinematic conditions can always be expressed by equations between 
the coordinate variables, preferably in such a way that the aforementioned 
desirable features are achieved. 


The materialization of such a kinematic equation between any two or 
three of the coordinates is the problem of the design engineer, who would 
have to decide whether to use linkages, gear wheels, cams, or hydraulic 
connections, and so on. This detail— designing problem is beyond the 
scope of this report. 
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Friction forces at the hinges, particularly at the lagging hinge £, 
have "been tried, for the purpose of damping excessive lagging oscillations, 
"but with the detrimental effect of producing high bending moments at the 
blade root. Such devices can also be readily calculated by equations ( 97 a), 
( 97 b), and ( 97 c). 

The four cases of kinematic and friction constraints previously 
enumerated will now be discussed in detail both for hovering (n_ * 0) 

and for small speed ratio (n e / 0). 


Case A: 0 (Pitch) Fixed, 3 (Flapping) and £ (Lagging) Free 

In this case, 0 «= 0 . The moments and M y1 will be zero, but 

not the twisting moment M^, which will be taken by the pitch-holding 

and pitch-changing mechanism or by counter flyweights. Two equations, 
namely (97b) and (97c) with the right-hand sides. equal to zero, must 
therefore be solved for two unknown variables 3 and X’ 

Neglecting I and I zg in comparison with I H (cf. equations (77)) 

and observing in accordance with expression (78) for B that in equa— 

tion ( 97 b) the Coriolis term and the aerodynamic damping term in X partly 
cancel each other, whereas inequation ( 97 c) the Coriolis term and the 
aerodynamic damping term in 3 can be combined as one term, the equa- 
tions of oscillation (97b) and (97c) became simplified to 


- Kg + 'f e 0 i 12 - 1 1 $■*£**&} + 



t I _ p £l i£_ 

h p a 315 


1 + t e iXf 




(98b) 


Equations (98a) and (98b) are a system of linear, homogeneous, differential 
equations with constant coefficients. The complete integral of these 
equations can, as is well known, be built up by particular integrals of 
the form 


3 - Fe 5 *, X “ 


(99) 


where F and D are real or complex constants (amplitudes) and Sr =r fit. 
The values of the real or complex constant (frequency) p must be deter- 
mined, as usual, from the .condition that equations (98a) and (98b) have 
solutions different from 3 = 0 and X * °. (This means that their 
solutions are also different from F = 0 and D = 0 .) These two equations 
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vill also determine the ratios D/F of the amplitudes of oscillation. 
Putting equations (99) into equations (98a) and (98b), noting that 

— a ft, and dividing through by ft^ yield the following homogeneous 
dt 

equations in F and P: 



As anticipated before, the determinant of the coefficients F and D 
must vanish. Hence p must be determined from: 



or 



The value of © c is given by equation ( 20 ). For example, if 

W s 4000 pounds, p = O.OO238 slug per cubic foot, E s 25 feet, 
ft = 20 radians per second, a_ » 0.020 (angle of attack about 0° for 
Clark Y airfoil). 
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ri e » 0.05, i x 25 * ^.17 feet, n * k blades, and * 0.2, (102) 


then the value of 6 C will be 


9 q m O.O3O6 


( 103 ) 


Taking the valve of — , moreover, from equation ( 77 ) s namely 

B 0.00386 

’ ^ 


gives 


/ -\ 2 

e ° W " 1-39 x 10 


~ - s - § e c » 1.52 x 10" 


Ms * 


2 S g 

T fiP 


0-618 x 10-5 


x 10 


(101O 


Consideration of the four quantities of equation (lOh) which appear in equa- 
tion ( 101 b) as negligible in comparison with unity is therefore Justified. 
The second and foprth terms of equation ( 10 h) represent the influence of 
the weight of a blade, and this influence will, in accordance with the 
foregoing considerations, be henceforth everywhere neglected. Equa- 
tion (101b) may hence be reduced to: 


p 4 + Tp 3 + (l + *%)p 2 + 2 T| e Tp + 2 q e ^l + 2 t) € 


( 105 ) 
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vhere 


T - l 32 
(1% 315 


(‘ 



1 


+ 1 %) 



Because of the assumption (se 


thickness t^ of a blade 




appendix G) of a fairly high relative 
it must be assumed in the calculations 


(c 

that the blade material is of vood. An average value of the ratio 
will then be (for p = O.O765 lb/cu ft, a » 30 lb/cu ft) 


p/a 


£ » 0.0025 

a 

(For a hollow, built-up blade design the average density ratio may be even 
somewhat smaller.) Hence, by using equation ( 77 )* the numerical value 
of T is found to be 


T =» O.595 


( 106 ) 


As can be seen from equation ( 101 a) , equation (105) can be readily solved 
by writing it in the form 



+ Tp + 1 + 2 q 


e) 


0 


(107) 


The solutions of equation (107) are: 


* 1,2 a±L \l*% 



These solutions, which are of the form 


P 



S 


i 


(108) 


(109) 


have, in accordance with equation (99) > the following physicsl signifi- 
cance. The logarithmic decrement log , defined as the logarithm 

A n+1 
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of the ratio of the amplitude of one cycle to the amplitude of the 
succeeding cycle, will "be 


log 


Ann 



(no) 


The natural (real) frequency of oscillation q in cycles per second will 
"be 


q = — S (111) 

2it 1 

where fl is in radians per second. 

Thus, the solution (108) with the numerical data of equations (102) 
shows that in case A there will he two natural frequencies: namely. 



1.005 cycles per second 




« 3 . 3h- cycles per second 


The oscillations corresponding to the higher frequency 
very highly damped, the logarithmic decrement being 



1 + ^ 


= 1.78 


q„ . will he 

3 , 1 - 


* 3 , 1 - 


to the rotational frequency 
(flapping) will, because 


Therefore, despite the proximity of 
— , the oscillations corresponding to 

2rt 

of the high damping, present no danger of resonance and will he quite 
Btable. (In order to get a more exact insight into the influence of 
res onanc e , it would he advisable to determine the ratio of thrust 
fluctuation to flapping amplitude. This determination requires a study 
of fluctuation of engine torque and speed. This question of resonance 


will not appear in the case (B) of appropriate kinematic constraint 
between flapping and lagging. ) On the other hand, the oscillations 
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corresponding to the low natural frequency q^ _ (lagging) will "be 
praotioally undamped. * 


The ratio of amplitudes can he obtained by putting equation (108) 

into equation (100a). Thus, corresponding to the complex frequency p , 

__ *2 



^c 

.a 

liil\S 


- 1.18 x lcT^l + O.lSli) 


(Equation (100b). because of the (nevertheless close) approximation. for p, 
would give zero.) Quite generally a complex value of F/D can be 
interpreted as follows: If 



and if, correspondingly. 


- « a 1 bi 


then, with the arbitrary amplitudes and Hg, either ? or p 
(say, X) will have the form 


X ■ e -Re ^jK L cos + Eg sin 

and, from the amplitude ratio F/D, P will be given by 


P 


a£(*) + b?Nr + 


2S. 


(112a) 


(112b) 


where ?(y) is the expression in brackets appearing in equation (112a). 
A complex value of a ratio F/D of amplitudes therefore indicates a 
difference in phase between flapping (p) and lagging (?) oscilla- 
tions, the real part giving the magnitude of the component of "ft in 
phase with X and the Imaginary part the magnitude of the component 
of p one-quarter of a period out of phase with X» 
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Because of the actual small absolute value of 


It is seen 


that the undamped oscillations, corresponding to the rather low natural 
frequency T i e > will occur practically only about the z^— axis; that 

is, only the lagging angle will oscillate. 

Following from the complex frequency p_ ,, the ratio D/F will, 

• 3 J 4 ** 

from equation (100b), be 


D > 

F/3A 


9 °fe p 2 3 A “ 315 



1 + 


t*l) (|° * °fe) S 3,l 


0.01005 ± 0.02951 


P‘ 


3,^ 


+ 2q, 


From the fairly small absolute value of the ratio (D/F) . it is seen 

3,4 ft / N 

that the highly damped oscillations of the natural frequency — (1 + r\ e ^ 
will occur practically only about the y^— axis; that is, practically 
only the flapping angle will oscillate. 


The results for both frequencies, implying that the lagging and 
flapping oscillations are practically independent of each other, show 
that the slow undamped lagging oscillation is very sensitive to dis- 
turbances and that the high restoring forces in flapping have no com- 
ponent which might oppose the lagging deviations. 


Case A^: 9 Fixed, 3 Free, 5 under Friction Constraint 

The unfavorable result concerning the lagging oscillation leads 
to the attempt to improve the stability by introducing a friction con- 
straint acting on the lagging angle by means of a moment K£ft at 

the root of the blade, where K denotes a constant of relative energy 
dissipation the value of which can be chosen according to the required 
degree of damping. This constraint could be accomplished, for instance, 
by a rod leading from the root of the blade to the hub and provided with 
a telescopic fluid brake. 

« 

The term -K£ must then be added to the left side of the dynamic 
equation of the lagging acceleration, that is, to equation (97c). Then 
by the s ame procedure as before the amplitude ratios (see equations (100a) 
end (100b)), the frequencies, and the damping terms (see equation ( 103 )) 
can be determined. It is considered sufficient here to calculate only 
the new complex frequency consisting in the damping and the real— 
frequency contributions. 
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Instead of equation (105), the following frequency equation appears: 
p^ + (T + E)p^ + ^1 + 4q^p 2 + 2t| e (T + K)p + Ep 


+ 




«= 0 


(113) 


In order to determine the effect of the constant K in general - 
terms on the damping and on the frequency, a simple approximate general 
solution was obtained by Newton’s method based on the assumption that E 
is small. The first approximation - here was taken as the exact solution 
to equation (113) when the term Ep is neglected. This first approxi- 
mation is then the same as the solution given by equation (108) with T 
now replaced by (T + K). The second approximation is then found to be 
the following: 


* 1,2 “ " 


P 3,4 - - 


where 


l»En 


+ ' 1 ~ if ~ * K jl ‘ 

[f - 4*, e (T + E)j 2 + 8rj e [e - 4n 0 (T + Ej| 2 + 8n 0 




/m ^ + K ) A + 2(1 + 2tJ e ) 

(T + E) _ K le 

2 


A 2 + 4(l + 2q e ) 


i I 


(l * 1.) 


> 


1 - 


A - T - E 


A 2 + 4 


( 1 + 2 0 


J 

(114a) 


A _ (T + E)3 

A = ■ ■ ■ — + 

4 


(T + E)(2 + 4tJq) + E 


For an example, the same numerical data as assumed in the preceding 
section (see- also appendix G) together with the value of E = 0.1 were 
introduced. Equation (ll4a) then gave the following results: 


p l,2 = “ 0,05 ±_ ^319i 
P 3 ^ = - 0-310 + 1.0351 


(ll4b) 
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The validity of the approximations given by equation (11 4a) was 
checked by putting the same numerical data into equation ( 113 ) and 
solving it exactly (by Ferrari* s method) to three significant figures. 
The results were: 


v 1)2 = - 0.0525 + 0.3211 

> 

p = - O .295 ± 0.9721 
3,4 


(114c) 


Comparison of equations (11 4b) and (114c) shows that with values of K 
not appreciably greater than 0.1 the comparatively simple formulas in 
equation (114a) are sufficiently exact for most practical purposes, so 
that it is unnecessary to formulate an exact, but more involved, general 
solution of the quartic equation ( 113 ) here. 

The natural frequencies and the logarithmic decrements corresponding 
to the solution (ll4c) are (see equations (110) and (111)): 


^ 1,2 


1.01 cycles per second 


q , = 3«30 cycles per second 



where, temporarily, denotes a real amplitude. 

The following conclusions for case k^ can be drawn from these results 

of the numerical example. It is possible to achieve fairly high damping 
of the lagging oscillations by introducing moderate fluid friction. 

This friction, moreover, will have little influence on the two values 

of the natural frequency and, on the high damping which is associated , 

with what originally were practically the flapping oscillations. It 
must be observed, however, that because of the low natural frequency 
(q-^ 0 ) corresponding to the lagging oscillations, the restoring force 

in these oscillations will also be low, and may, in fact, not be 
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sufficient to return the blade to its normal, position. The damping intro- 
duced by the fluid friction would in that case be of no avail. 


Case B: 9 Fixed, 0 and £ under Kinematic Constraint 

It has been seen (case A) that with the pitch angle 9 fixed and 
the flapping and lagging angles P and £ free, there will be the 
disadvantage of an absence of damping for lagging oscillations. The 
possibility of overcoming this disadvantage by introducing kinematic 
constraints between the variables will now be investigated. 

First an appropriate constraint between flapping and lagging will 
be introduced. As previously explained, this constraint should be such 
that the conditions of steady flight are not violated; that is, the 
constraint, as represented by an equation, should satisfy the condition 
that when £ =■ f^,, then 0 = 0 C , or when £ = £ — £ c = 0, then 

f s P - p. * 0. This condition will not be violated by a constraint of 

G 

the form 


£ - *3 


(115) 


where k is a constant to be chosen in accordance with requirements of 
stability and of avoidance of resonance. In order to achieve materially 
such a constraint, a preadjustment of the angle 0 or the angle £ is 
necessary. For example, the lagging angle £ may be preadjusted to a 
value where £^ is the value of £ when P => 0. From equation (115), 

this value is ^ *= £ c — K P c . 

As in case A, the solution for 0 and £ will be of the exponential 
form given by equations (99) « When equations (99) are then put into equa- 
tion (115), the ratio of the amplitudes is obviously 


D 

F 


K 


(116) 


Condition (115) can mathematically be taken into account by means of a 
Lagrange multiplier and a Lagrange function 0, where, according 

to equation (115), 


f = £ — K0 ss 0 


(117) 
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The equations for 33 a nd F corresponding to equations (100a) and (100b) 
then become: 


+ Dpe c I + Lm 


j ( p2 + 1 + S + 

,?9 ° - p ^ IX + + ap )l ■ + ^ 


+ L M a 0 ( 118 a) 


S3 


0 (ll8b) 


From equation (117)# 


^ «= — k, and ^ » 1 

S3 


(119) 


Putting equations (116) and (119) into equations (ll8a) and (118b), and 
eliminating 1^ from these equations gives 


F < 


“ rj2 - - 1 1 - P — “ 


“ %( ] 


p* + 1 + 


fl 2 315 


- I s ! + I 1 ®) 


+ p K 0 O I 


+ * 


“ P JjS + i^(f S = + “t)l " + **') | ’ ° 


( 120 ) 


For equation (120) to have a solution other than F =» 0, p must have 
the values satisfying the equation 


— p 2 (l + k 2 — 2 k 9 ]— pa* 


- |i + 2n e ^i + 


0 ( 121 ) 


where 

a » 32_ 

" n 2 ?- 315 




f 1 + K + l>°) 


Hence > neglecting the quantities 2 3. — in comparison with (l + 

r 

and also neglecting a* 2 in comparison with 4^1 + k 2 ) 1 + 2q e ^l + k 2 ) 
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P ~ 






( 122 ) 


It may Ids observed that In this case of geometric constraints, the drag 
has a slight, hut noticeable, effect on the damping; whereas in the 
case of free oscillations, the Influence of the drag on the damping 
was negligible. For ordinary values of k, the damping will be quite 
large, and the stability therefore very great. For example, if 


K s 1 


then 


p --r 1 /!* 1 < 123) 

By use of the same numerical values as before (equations (102) and (103)) 
there is obtained 


p =« - 0.158 ± 0.775i 


(124) 


PO 

which gives a natural frequency of q =» •=— O .775 =■ 2.47 cycles per second- 

2 it 

and a logarithmic decrement (see equation (110)) of 2 jt - * ■•- ? ■ * 1.282. 

0.775 

From the point of view of stability this case therefore appears quite 
satisfactory because there is considerable damping and there is little 
danger of resonance from the drive, that is, from such probable dis- 


turbing frequencies as n 


( 20 
\ 2 « 


a 3.19 cps j or 20 . 


Case Cj 9 , 3 , and £ under Kinematic Constraint 


In the cases thus far discussed, the pitch angle 9 has been 
assumed fixed. The effect of allowing a certain freedom of pitch 
change 9 in accordance with geometric conditions (constraints) among 
the angles 6 , 3, and £ will now be considered. For this purpose, 

constraints similar to that used in case A (equation* (115) ) will be 
assumed, namely. 
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P 

I 



(125) 


where X and p are constants which will he selected in accord with 
stability requirements as will he discussed. , 

Solutions to equations (97a) to (97c) will have the same form as 
given hy equations (99), namely. 


W « Be p *, ? = Fe^, £ « De 5 * 


(126) 


From equations (125) and (126) it follows that 


F a XB 
D a pB 


(127) 


By putting equation (126) into equations (97a), (97b), and (97c), equa- 
tions of the following form are obtained: 


BP a> + re er + Dp 2 d 


BP- V + FP_- + DP 


M yl e 






> 




3h 3 f " 3d ~ 


(128) 


where P indicates polynomials of second and lower degree in p. As in 
case B, the constraint conditions (equations ( 125 ) ) can he taken into 
account hy means of Lagrange multipliers and and kinematic 

conditions 0^ and (%, where 


0 1 a p - xe a 0 

02 a 5 - 1^*0 
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Eliminating and from the three expressions of equation (12<?) 

results in the following single equation: - 

B j^lb + XP lf + ^ P ld + + ^ P 2f + + "( P 3 h + ^ P 3f + ^ P 3d)j * ° 

The value of p is then determined by setting the factor of B In 
equation (130a) equal to zero. Equation (130a) can hence be written as: 

P lb + X ( P W + P 21>) * l '( P ld * P 3l) + x2p 2f * “Si + X K P 3f + P 2d) - 0 < 130b) 

When the expressions obtained for P (see appendix H) are substituted 
In accordance with the definitions (equations (128)) and the expression (78) 
for as well as the geometric properties (equations (77)) of the blade 

body, are used, equation (130b) becomes the following quadratic in p: 
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With, the numerical values consistently used in this report, equation (131) 
reduces to (see appendix E): 

2 

p + 2 a.jP + a Q » 0 (132) 


where 


» (A 2 + u 2 )"" 1 (0.299\ 2 - 0.0068X + 0 , 00126\|j, - 0.0082 m. - 0.00084) 

a 0 « (A 2 .+ U 2 )"’ 1 (1.1X 2 - 0.560*. + 0.1m 2 + 0.00656) 

The solution to equation (132) is 


P 


a l 





(133) 
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For stability it Is necessary and sufficient that the roots of equa- 
tion (132) “be either real and negative or complex with the real part 
negative. These conditions will he satisfied if and only if a, and a 

o 

are both positive (or zero). The values of X and u can now be bo 
chosen that these conditions will be satisfied. For example, if 


\ = 0 . 7 , 4=1 


then the complex relative frequency p will be 

p = - 0.0882 ± i x 0.403 (134) 

0 0882 

giving a logarithmic decrement of 2 jt — * 1.370, and, with 

0.403 

n = 20 radians per second, a natural absolute frequency of 
20 

— 0.403 « I.287 cycles per second. From the point of view of stability, 
2?r 

case C, as well as case B, is therefore satisfactory for hovering. 

OSCILLATIONS OF BLADE SYSTEM IN LOW— SPEED TRAVELING 


The same assumptions in regard to the (first order small) value of 
the speed ratio u and the neglect of third-order terms, as made for 
the steady state in the section STEADY STATE IN HOVERING AND IN LOW-SPEED 
TRAVELING, will also be made in this section for the oscillations in 
traveling. 


General Explicit Equations 

As in hovering, the oscillation equations ( 74 ) to (76) must be 
solved, except that now the expressions for 0 Q , p Q , and £ 0 will be 

not only different from the hovering values 0 C , 0 C , and C c , but also 

variable, and in addition the equations will contain more terms than in 
hovering. 

The expressions for 0 O , 3 o , and are given by equations ( 84 ) 

with- (96), (85), and ( 87 ). By using equations (35) to (37) and ( 49 ) 
to (?1) for the quasi-elastic and the damping momenta and rejecting all 
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terms smaller than, the second order, the dynamic equations (74) to (76) 

( _ Jtp k p \ 

E = — B as in the foregoing section J become : 


“ "^dzs + 



Ur) h( 2 + + 7T1 * + 8111 ^ 


I (f — 20 - Icy I — ) + - f ~(l + 2T| + an e sin ♦) 

ft l B c l^ac j j® / ft R \ e ® ' 

E(l + s?) c, - B(l + s,) c, 

_■ '- ■ < -if— ftX„ + — -i 


12 E 


E(l + Si ) c ± ~xl 

— 3 — — - fu cos * = 

6 E e ¥ <X 


3 B f ^e sin + 


M_ 


(135a) 


*A. “P(^ + V) + ^0^2 + + 

(l + |s^E(l + 3d e Bin *) - pjft 2 ^ + Se) - ~_(l + |s^ 


+ 0 32 


Eu e cos i(f 


(135b) 


- *«’ - * 9 As " + i M ) - 1 


+ — Ekt +-0]0- £ft 2 Se = 

315 \ P 5 V or 


2P | 32_(l + is )(| 0 - a 

c ^ ft U5\ 2 V\5 c P ; 


(135c) 


When equations (135a) to (135°) sre compared with the corresponding set 
(equations (97a) to (97c)) for hovering, it is seen that both sets are 
the same except for four additional terms in equation (135a) for 

and three additional terms in equation (I35t) for M^. Each of these 

additional terms has the periodic coefficient sin ijr or cos » nt) 
but of relatively small magnitude. It will therefore not be necessary 
in this case to apply the theory of Mathieu functions, but it will be 
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sufficient to use a method of successive approximation starting from the 
hovering state as a first approximation. This method may be expressed 
by the propositions 


9 = 0 1 + 7 - % + 4 e 3 2 , r - fi + ^ e r 2 


The validity of this method of solution vill be checked by the 
results obtained by its application, for it is necessary that these 
results remain within the order of magnitude assumed in advance; that 
is, in this case the terms proportional to u e ^e.g., u Q 9 2 j must not 

be larger than second order small, so as to be of a higher order Bma.ll 
than the corresponding hovering solutions ^e.g., 9^. 

Case A: 9 Guided, 9 = 9 Q (ijr); 3 and t. Free 

This case is the same as case A for hovering, except that now instead 
of keeping 9 fixed at the steady-state value 9 C for hovering, the 
blade angle is guided so that at all times 9 « 9 Q , where 9 Q is a 

function of ^(^ftt) and therefore of the time and is given by equa- 
tions (84) and (96). The deviation of 9 from the Bteady-state values 
will therefore he zero, so that, as for case A in hovering. 


9 = 0 

and only two equations, namely equations (135b) and (135c), need be 
considered to determine the natural frequency and the damping. 

Equations (135b) and (135c) can be solved as indicated above by 
writing the solution in the form 


6=0, 3 = + h e P 2 > r “ {>2 U36) 


where P 1 and are the solutions for hovering, which have already 

been obtained in the section OSCIUATIONS OF BLADE SYSTEM IN HOVERING. 

Inasmuch as equations (135b) and (135°) are a set of differential 
equations of the second order in the two unknowns F and , £*, the 
complete solution to these equations must contain exactly four 
arbitrary constants, which must satisfy any given initial conditions 
of position and velocity. Moreover, any solution containing four 
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arbitrary constants and. satisfying equations (135b) end (135c) will "be 
the com ple te solution of those two equations (if it does not violate 
the condition that the results remain within the order of magnitude 
assumed in advance). Now, from equations (99) and ( 108 ), 

h - i v Pn ’ 1 '- h - ^x• sEn ’ ,, (K7) 

n=l n=l 


where F n is an arbitrary constant, whereas F n * is a constant depending 

on F . Therefore, from equation (136), it follows that, to obtain a 
n 

complete solution for p and £ in traveling, it is sufficient to 
obtain only a particular integral of the differential equations (138b) 
and (138c) in the unknowns and. £ . In obtaining a particular 

integral it will be sufficiently exact to deter m i ne P 2 and £ 2 to 

only first— order small quantities, because in the final solution they 
must be multiplied by the first-order small quantity n e and will 

therefore yield only second-order small terms. 


Thus, putting equations ( 136 ) into equations (135b) and (135 °) , 
dividing through by n e I H , using the expression (equation ( 78 )) for P c , 

and rejecting all terms smaller than the first order result in the 
following equations: 


- "s " + h * K) - V m f 1 + I s 


^1 

+ — sin if 

n 


— p^ cos + 


to * 0 


(138b) 

(138c) 


A particular integral of equation (138c) is obviously 


£2 * 0 


( 139 ) 


A particular integral of equation (138b) can be obtained by first 
observing the following relations: 
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s — = fi i- 
dt <i ♦ 


cos = 


,*♦ + e" 1 * k 


(140) 


sin t - - — (e 1 ^ - e" 1 *) 


By putting equations (137) and (140) into equation ( 138 b), the equation 
for 3 0 becomes: 


P 2 " + T3 2 » + 3 2 



where 

• -L. 


T s 


E 


0^ isC 1 + h * K) 


*= 0 . 595 . (Cf. equation ( 105 ).) 


A particular integral of equation (l4l) will be the sum of four pairs 
of terms, each pair corresponding to a given p n « Thus, for each p n . 


W -A n .M* + A , e M* 

' ^n n 


(142) 


The ratios of the constants A n and A n * to F q are readily obtained 

by putting equation (142) into equation (l4l). Thus, by writing 
(cf. equation ( 109 )) 


S \n + 1S ln 


the value of A n is seen to be given by 
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4 A n 


1 + 




^-l"e)* n + *)* + T ( P * + *) + 1 


C 1 - S ln) 


+ IB. 


en 


Sen 2 - ( s in ♦ 1 ) 2 + ™en + 1 * 1 ( S ln + X ) (“en * T ) 


Eat ionali zing the denominator yields 

^ - fi^K 

where 

C « E 


rv v ~ 1 



\ 

(l - Sln )G + R en D 

+ i 

3»° - ( 

1 - SmjD 


2 2 
C + 


(143a) 


+ ™en + 1 


Similarly, 


V = iK 


- (v + 2 ) 

D ■ ( S ln + ^(fen + T ) 

K 1 t s to)°' - ~ 1 _B.„C* + C 1 * 3 to) D l 

0*2 + d*2 


(143b) 


where 



Taking the values of R en and S ^ n from the hovering solution (equa- 
tions ( 108 )) and the definition (equation ( 109 )), and observing that t he 
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numerical calculations are simplified "by the fact that, for n => 1 , 2 , 
\n 

of A n /F n and V/V 


E = 0 and, for n =» 3»4, D = D* = 0 , yields, for the numerical values 


Sy^i ■ - a - A i*/ F i = c + 4i* 


■ c — di, A q yF^ ■ — a + hi 


2/2 


A y ^3 53 ® + fl > A 3 y ^3 ■ S + hi, Ay^ - g - hi, - e - f i 


(144) 


where from equations (143a) and (143b) 


a » 0.182, h = O.I 97 , c « 0.641, d * 0.486 


e = 0.0064, f * 0.0348, g « 0 . 947 , h > 0.126 


The numerical values of equation (144-) show that Pg actually is a 
first-order small quantity ^that is, of the same order of magnitude as p-jJ 

and that therefore the method of solution used here 1 b valid, inasmuch as 
P 2 always appears in the form a Q F 2 . 

The foregoing results can he mathematically interpreted for each 
pair of conjugate values of p Q as follows: 

With 


P 1 " ” E el + s U i » p 2 * “ K el “ S u 1 


for example, and with a, h, _c, and d defined as in equations (144), 
it can he easily shown that f p. ) is of the form 

V -‘-'1,2 






(145a) 


where 



NA.CA TN No. 1^30 


69 


"l.aCvO * ^ ° oe (V) + H 2 sln 



n-nri and are constants. 

_ It can be. readily proved, moreover, that the additional solution 
(P 2 )l 2 will then have the form 

(Pg) = e_Bel1,f ” a ^l,2 (pil + lj * ) + c ^l,2( S il “ X * *) 

1,2 


bp. 


i,2 °n 


+ 1, Iff + 


2 ( s il 


+ 1 


S + d Plj2 fs ±1 - 1 ,+ + - 


(sn- 1 ) 


(Cf. equations (112a) and (112b).) 

Comparison of equation (l4-5b) with equation (145a) shows that a 
complex value for the ratio of amplitudes indicates a difference in _ 
phase between hovering vibrations P, and the additional vibrations P Q . 

1 / \ jj / \ n * 

In fact, for each of the frequencies fs^ + l) — and — l) — , 

(^ 2 ) consists in two components one— quarter of a period out of phase 
with each other. 


The physical significance .of the results can be stated as follows: 
In regard to the frequencies, equations (142) and (108) show that one 
new frequency is added which is practically double the frequency of 
rotation 2. Three other frequency roots appear which, howe ver , are not 
significantly different from the frequencies of the roots \f&% 7^ an<i 

^1 + of hovering. 

The additional amplitudes given by equation (144) in terms of the 
amplitudes in hovering, and depending on initial disturbances, are small 
in comparison with the hovering amplitudes, as long as the value of u e 

is smal l. This result, in fact, is the proof that the method of integra- 
tion is consistent. 


The logarithmic decrements remain practically unchanged in the 
transition to traveling, and the solution corresponding to the frequency 



which has no damping in havering, therefore is still very 


sensitive against disturbances in the transition to traveling. 


(145b) 
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Case B: 0 Guided, 3 and £ under Kinematic Constraint 

It was seen that in case A of traveling there is the same danger of 
lack of stability due to the absence of damping as in case A of hovering. 

___ Case_B of hovering vith-an appropriate kinematic condition between 
3 and £ made it possible to obtain better damping and also frequency 
values of no resonance danger. Whether a device of the same kind will 
serve the same purpose in traveling will now be determined. 

For case B of u Q 0 with its kinematic constraint between 3 and 

£ it will also be advisable to avoid mutual bending moments in the con- 
straint mechanism in the steady state of motion by a pre adjustment 

3_„ — (3)„ _ or £ = (£) . (During the oscillation, however, such 

pr £«=0 pr 3=0 

mutual moments, though small, can again not be avoided. ) Such a pre- 
adjustment for the case of traveling must, however, be periodic, as can 
be seen by the following consideration. The kinematic condition may 
again be expressed by 


T * *3 or £ — £ q « k(3 — 3 0 ) (146) 

where k is a constant and £ 0 (i|r) and 3 0 (>|f) are periodic functions 

of the angle of position +, given by equations ( 85 ) and ( 87 ). From 
equation (146) it follows, then, that the preadjustment must be 

3p r » 3 q ~ (equivalently) £ pp = £ q - k 3 q . Some means, for 

instance, a cam plate, will be required to enforce such a periodic condition. 

The constraint (14-6) can be treated, as in the case of hovering, 
by means of a Lagrange multiplier L . Thus, equations ( 135 b) and ( 135 c) 
become: 


+ -HA- sin i 


■) 


- + v v - 1 1 

+ ^) - ^9 »»■ *J “ v 0 


- 3 


3 9 C I 


- I I H - 3 23 c «E h - | ||^ (l + e c ~ a p J 


- la% + 


0 


( 147 c) 
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Eliminating 1^ from equations (147b) and (147c) and thus obtaining a 

single equation, expressing f ty P according to equation (146), 
setting 

P - ?! + d e P 2 (148) 


vhere 3^ is the solution for P in hovering (fx e = 0), rejecting terms 
smaller than the second order, and finally dividing through by 4 e T H 

result in the following differential equation for P 2 : 


- p Q (l + k 2 ) - Tpo - P 0 fi 2 



sin Mr — 3^ 


cos 


♦) 


(149) 


where T is as defined in equation (105). It will be observed that 
equation (149) Is similar to. equation ( 138 b) of case A except for the 
coefficient (1 + k^) of Kj in equation (149) and for the fact 
that p-]_ must now be taken from case B of hovering instead of from 
case A. Thus, from equations (99) a 11 *! (123) j 


P, 


F n e 


n t 
F n 


( 150 ) 


For a complete solution for (3 it is sufficient, as explained in 
the previous case, to obtain only a particular integral of equation (149). 
This particular integral can be obtained in the same manner as shown in 
case A for equation ( 138 b). Thus, 


AtS 


+ „( P2+1 ) 


+ A g e 


+ A 2 *e 


( P 2 _1 ) 




(151) 


where A must be determined in terms of F„ by substitution for (3 

n v 2 

into equation (149). By use of the expressions for sin cos t, 
and (equations (140) and ( 150 )), equation (149) can be written in 

the form (cf. equation (l4l)): 
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When equation ( 151 ) is put into equation (152), is is found, for 

. ft = 1 

(which gave satisfactory results in regard to stability in hovering), that : 


where 


V 


■ T- T 




(±- 

3_V i 1 " 

S ln)° + E en B ] 

1 

1 Pm 

•H 

+ 

1 - s i>_ 

V 


0 2 

+ D 2 



0 - 2V 2 + 

™<m + 1 “ 

3 In + X ) 2 



E * (“an + 

T )( 3 in + *) 



3 ^ 

1 , £ i+s ^ 

1°' - H en B '] - 

‘[v 0 * ♦ k 

♦ s in) B : 

2%, 

rn 

c* 2 + 

D* 2 



(153a) 


(153b) 


where 


°* - + «*«, +- 1 - S ( S ln - 1 ) 2 

D- . (l® 0n + l)(s to - l) 

With the values of R Qn and given by the solution (l 2 h-) for 

havering, equations (153a) and (153b) yield the following results: 


A = F-^C— a + bi), A 1 * » F^c + di) 
A 2 = F 2 (c - di), A * » P 2 (- a - bi) 


(154) 


where (from equations (153a) and (153b)) 
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a - O.OI75, b - 0 . 012 l )-5 

O «* 0,860, 4 a O.O678 


These results can "be physically interpreted in a manner similar to 
that in case A. In the transition to traveling, two new natural fre- 
quencies are added to the frequency O.775 for hovering. These frequencies 

are (1 + O.775} — and (l - O.775) — or 1.775 — and 0.225 — • 

2 « 2k 2k 

Neither of these frequencies appears to present any particular danger 
of resonance, which might otherwise "be caused hy exciting disturbances 
having a frequency connected with the frequency of rotation, that is, 

a multiple of The damping decrement, moreover, remains practically 

27C 

the sane in low— speed traveling as in havering (log decrement * 1.222), 
and this means that the rotor system will remain quite stable in the 
transition to traveling. 

From the fact that the values a, b, c, and d in equations ( 153 a) 

a n d (153b) must be applied with the small factor u in accordance 

© 

with equation ( 148 ), it is seen that the additional amplitudes in the 
transition to traveling are again small in comparison with the amplitudes 
in the state of hovering. Thus, it appears that case B (of constraint 3 , £ ) 

remains satisfactory in the transition from hovering to traveling. 


Case C: 0 , 3 , and £ under Kinematic Constraint 

This case has (as in hovering) the practical advantage over the 
others in that the blade a n gle need not be guided but will automatically 
adjust itself to the proper value to support the weight of the heli- 
copter. The s a me constraints for 3 and £ as in hovering will be 
assumed here, namely 

3 » xe 

r - 4-is 

where X and u are constants. 

These constraint conditions can be realized practically, for a 
correct pitch-angle function 0 Q , by means of a variable pre adjustment 

of the angles 3 and £ . This preadjustment is given by 
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(156) 


These pread,}ustment values can he found explicitly hy substituting the 
expressions for Q q , 0 q , and £ Q given by equations ( 84 ) with (9 6), (85), 
and (87) , respectively, into equation (156) . They will appear in the form 

Ppr “ f 1 + V e (t 2 + f 3 Bin * + f ^ 00s 

5 pr s + ^ e ( d 2 + d 3 8in t + d 4 oos ^ 

where f and d are constants found from equations ( 84 ). (96) (85) 

and (87). 

Thus (3 and £ will be fairly simple functions of the angular 
s? L pr 

position i|r of a blade, and can be materialized by, for example, a cam 
plate . 

The stability of the rotor system with the constraints ( 155 ) can he 
considered in the same manner as case C in hovering. Thus, using 

Lagrange multipliers and L^ and denoting by D 1 (N', 0, f, if), 

D g , and the left-hand sides of equations (135a), (135b), and (135c), 

respectively, yields, for these equations: 

aj t 

V 0 , T, t) + L^i -i + -a „ ° (157a) 

, 30- 300 

D 2 ( 0 » +) + — i + L^ — - ° ( 157 b) 

OP op 

30. 50 

D 3^ * ^ « + \l + ^ ^ ■ 0 (157c) 

where 0 1 and 0 g are the functions given by equations (155). Substi- 
tuting the values of the derivatives of 0 1 and 0 g , solving for 

and L^g by means of equations ( 157 b) and ( 157 c), and then substituting 
into equation ( 157 a) yields the following Blngle equation: 
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D 1 + XD 2 + hD 3 - 0 (I58) 

By taking the expressions for and from the left-hand 

sides of equations (135a) to (135c), eliminating £ and tf from 
equation (158) by means of equations (155) j setting 


0 


e i * ‘V’a 


(159) 


where 0=0^ is the solution of equation (158) in hovering, rejecting 

all terms in equation (158) smaller than the second order (remembering 
that u 0 is assumed first order small) , and finally dividing through 

by u 0 T h the following differential equation is obtained for 0 2 : 


Va + Vi fl + V/ “ =r \ + ¥0) sin * + W cos * 


E (1 + s. 


(160a) 
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The expression for 0^ is 


' 9 1 = 


B e 
n 


P n * 


(161) 


where p n is given by equation (133) and B n is an arbitrary constant. 

As already explained in cases A and B, it is sufficient in solving 
equation (l60a) to obtain only a particular integral of this differ- 
ential equation, and this integral can be obtained in the same manner 
as in cases A and B. Thus, by using equation (l6l), noting that t = ^t, 
and also noting the expressions (140) for sin y and cos t aB expo- 
nential functions, equation (l60a) may be written in the form 

(where * 

(v 1 )* 


a 2 02 



(f A - *0) 


+ b * 
o 


e 


+ b 


>) 


b * 

o 


(PA)' 


(162) 


As in cases A and B, set 


nil — 


(163) 


where A n end A n * must be determined in terns of B n by substitution 
of equation (163) into equation (162). Thus, 


A n = 


fl_ B 

2r n 


V + 1 (Bn' b l ~ b o) 


2a*T 


■H 


(Pn + O 2 + ( p n + *) a l + 


(164) 


By putting P n = B en + iS^, the following expression for A n can be 
readily derived from equation (l64) : 
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For the special case that was treated in case C (hovering), which 
gave satisfactory re stilts in regard to stability, the values of X 
and M. were taken to be 


X = 0.7, m, = 1 

With these values of X and u and with the numerical data (summarized 
in appendix G) which have been consistently used in-this text, the values 
of the constants appearing in equations ( 165 a) and ( 165 b) are, in accordance 
with equations (l 60 b), found to be: 


a 2 = - 1 .^ 9 , a 1 = — 0 . 280 , a Q = — 0.108 

E 1 

b, = — 4.55, b = — 13.22, b * = - 4.63, ~~ - 0.0417 

From equation (134), p^ is given by: 

p n = - 0.0882 ± 0.4031 

Hence R el = R e2 = - 0.0882, S 1]L = 0.403, and S ■ - 0.403. 

Substitution into equations ( 165 a) and ( 165 b) leads then to the following 
results: 

A 1 = B l^“ a + ti), A 1 * 

A 2 = B 2 (- c - di), A 2 » 

where 

a = 0 . 102 , 
c = 0 . 308 , 


B 1 (- c + di) 


B 2 (- a - bi) 


( 166 ) 




b = 0.0372 

d « O.O 695 
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As explained in cases A and B, the fact that the amplitude factors of S^e 

are of the same order of magnitude as the amplitudes in hovering shows 
that the method employed here of obtaining a solution in traveling is 
valid. 

The results can he physically interpreted again in a manner quite 
analogous to that in case A and in case B. Two new natural frequencies 

appear* in addition to the hovering frequency O.I1O3 - 2 -, namely, 

(1 + O.I1O3) — and (l — O.I1O3) — , neither of which appears to 
2 rt 2 jt 

present any particular danger of resonance. The damping decrement 
remains essentially the same as in hovering, and this shows that the 
system in case C will be equally stable in low— speed travel and in 
hovering. The natural inodes in this case can be expressed analogously 
to the relations (equations ( 145 a) and ( 145 b)) of case A (traveling) . 


CONCLUSIONS 


Following the ge one try, the statics, and the dynamics of the motion 
of a hinged blade system, the parameters of pitch angle 0, flapping 
angle 0, and lagging angle £ in the hovering and the traveling (that 
is, forward, backward, and sideways) steady state of flight were deter- 
mined. The geometric part of this problem, particularly for the case 
of traveling, consisted in the determination of the angle between the 
direction of flight and the zero— lift lines of the blade sections, which 
are rotating on a conical surface the axis of which is tilted toward the 
direction of flight. This involved the determination of the components 
of the relative inflow velocity both in the planes of the cross sections 
and in the direction of the blade axis. 

The influence of the induced inflow on the total inflow velocity was 
calculated only in regard to the direction of this total velocity, and 
the always very small induced change of magnitude of this resultant velocity 
was neglected. The three-dimensional Kutta-Joukowski theorem was then 
applied to the calculation of the lift— force vectors and their periodic 
deviations from the- planes of the cross sections. In this way the velocity 
components both along a blade axis and in a plane perpendicular to the axis 
were taken into account in thp vector product T x T, whereas the value of 
the circulation r was determined by the transverse velocity component 
only. In all previous publications the radial (blade-axis) component of the 
velocity had been neglected. 

The pitch angle 0 was expressed for hovering as a constant 
depending on the total weight of the helicopter and for traveling as 
a fractional function in terms of the first and secqnd powers of the 
speed ratio p. e , of the sines and cosines of the circumferential 
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angles and 2^ of blade position, and of the acceleration of the 
flapping angle 0. 

Although- the effect of the induced downwash on the lift forces has 
not been treated explicitly, this effect can, in accordance with the 
assumptions made in the present analysis, be considered to be contained 
implicitly. The induced angle a^, in fact, as given by equation (15) 

can be written in the form a« = — a. Therefore, if substitution of 

1 5 

(a — a^) for a is made in equation (7)* giving the magnitude T of 

the circulation, it follows that the effect of the induced angle on the 
lift loads will simply diminish these loads by a constant factor, 4/5. 

The numerical results given in this analysis will then actually remain 
unchanged if one assumes the gross weight V of the helicopter to be 
four-fifths of the value originally assumed. 

In the derivation of the steady— state values of the flapping and 
lagging angles based on the equilibrium of .moments about the hinges, 
the inertia moments of the angle accelerations 0, and £’ also had 

to be taken into account for the case / 0. The damping moments pro— 

, • • 

portioned to 9 , 0, and £ were also considered. This required the 

integration of differential equations in order, to determine in the section 
STEADY STATE IN HOVERING AND IN LOW-SPEED TRAVELING the steady-state 
values of the blade-position angles. 

The forces and moments due to small oscillatory displacements, 
velocities, and accelerations, necessary for the analysis of small 
oscillations about a state of steady motion, were determined in the 
section INERTIA FORCES AND MOMENTS AND EQUATIONS OF OSCILLATION. The 
inertial moments especially were expressed by means of the moment of 
momentum vector, and the Coriolis forces were obtained by the use of 
a rotating reference system. In this way, the complete system of the 
equations of small oscillations about a state of steady motion was 
established (equations (74), (75), and (78)). In the hovering state 
this system of differential equations has constant coefficients but— 
in the traveling state the coefficients have periodic additional terms. 

The integration was performed first in general terms, with results 
for frequencies, logarithmic decrements, and amplitude ratios given by 
simple formulas. These results were then applied for the following set- 
of plausible numerical design data. Four different cases in hovering and 
three corresponding cases in traveling have been discussed: 

Case A. Pitch angle 9 fixed, flapping angle 0 and lagging angle £ 
free. 

Case Pitch angle 9 fixed, flapping angle 0 free, lagging angle £ 

constrained by fluid friction (dashpot). 
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Case B. Pitch angle 9 fixed, 3 and £ connected "by a frictionless 
kinematic constraint. 

Case C. 0, 3, and £ externally free but internally connected by two 

(frictionless) kinematic conditions. 

The numerical data for these cases were assumed as follows: 

Total weight W = lvOOO pounds 

Tip radius R = 25 feet 

Rotational Speed SI = 20 radians per second = — =(3.19 cycles per second^ 

A - s \l/2 25 2rt ^ 

Chord c as cj — — \ r- = h2 


V- 


c i - 


feet 


Number of blades n 


Inner cross section of blade at -=■ = s. = 0.2 

R 1 

0 

Hinge eccentricity — = Ti e = 0.05 

*ku -I 

Thickness ratio of cross section of blade — ~ = Constant = — 

o 8 

Average density ratio of air to blade material — = 0.0025 

a 

Parasite drag angle ctp = 0.02 

In case A It was found that the oscillatory motion of the rotor 
system can be considered as consisting approximately of oscillations of 
only the flapping angle 3 and of independent oscillations of only the 
lagging angle £. The natural frequency of the flapping oscillations 

is = — (i + = 3.3^ cycles per second. Although this frequency 

is quite close to the rotational frequency — = 3.19 cycles per second, 
• 

there will be little danger of resonance because of the high logarithmic 
decrement, namely JtT/(l + q e ) = 1.78, associated with the flapping 

oscillation. 1 The lagging oscillations wili have the low natural 
frequency q, _ = — \/2r\ = 1.005 cycles per second but will be 

± » ei 2:t v e 

practically undamped, and therefore sensitive to disturbances. Inasmuch 
as the flapping and lagging oscillations are practically Independent of 


^or definition of T, see SYMBOLS. 
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each other. It follows as plausible that any phase difference between 
lagging and flapping which might arise on account of a separation of 
the flapping hinge from the lagging hinge would be quite small. It 
can, moreover, be observed from the formulas that both the natural fre- 
quencies and the logarithmic decrement are only slightly affected by 
the location t] q of the hinge, especially in flapping. 


In case A^, the stability of the lagging oscillations of case A 

is very much improved. by the introduction of fluid friction, producing 
a damping moment klgfl£ at the root of the blade; k is a constant 

of relative energy dissipation which may be adjusted to suit require- 
ments of operation. For the case of k = 0.1 it was found that the 
new natural frequencies will be practically the same as those in case A, 
whereas the new logarithmic decrement corresponding to the flapping 
oscillations of case A remains practically unchanged. The logarithmic 
decrement corresponding to the lower natural frequency, however, is now 
no longer zero, but fairly high. The numerical results for the natural 
frequencies and logarithmic decrements were: 


q, = 1.01 cycles per second, 

q^ h ~ c y cles p® r second. 



1.029 


1.91 


For any other values of k and any other numerical data, the results 
can be obtained by determining the complex frequencies p from either 
the biquadratic equation ( 113 ) or the approximate general solution (lit). 
The logarithmic decrements and natural frequencies can then be determined 
directly from equations (IO 9 ), (110), and (ill). 

Because of the low natural frequency, with the consequently small 
restoring forces, of- the independent lagging oscillations, the friction 
damping may prove insufficient for stability. As will be seen in the 
following cases, however, the damping can be successfully enforced with 
an appropriate kinematic constraint between, for example, lagging and 
flapping. 

Case B was worked out in detail for a kinematic constraint 
(geometric condition) of the form 
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where k is a constant for bttw.I1 oscillations. In cases B and C the 
method of Lagrange multipliers has "been used to satisfy such geometric 
conditions. These multipliers also have a physical significance, for 
they give the forces acting in the constraint connections. 

With such a constraint there will he only one natural frequency. 
(See equation ( 123 ).) For k = 1 and for the foregoing numerical data, 
the natural frequency of oscillation was found to he 2.k7 cycles per 
second, with a logarithmic decrement of 1.282. From the point of view 
of stability and avoidance of resonance this case appears quite satis- 
factory. For any other data, hut for the same form of constraint condi- 
tion, the frequency and the damping can he determined either by solving 
the quadratic equation (121) or by substituting in the approximate 
general solution to this equation (equation (122) ) . 

It may he noted that in cases A, A^, and B the aerodynamic loads 

had practically no influence on the natural frequencies of oscillation. 
As may have been expected, moreover, the natural frequencies are only 
little affected by the damping terms. 

Case C has the advantage that here the pitch angle 3 is auto- 
matically controlled. The two constraint conditions were assumed to 
be of the form 

p = X0, £ = M-0 

where A and p, are constants. This condition could be realized by a 

pre adjustment of flapping and lagging angles to the following values: 

B = 8 — A6 , £ =t— p0_. The method of Lagrange multipliers led 

pr c c 3 pr c c 

to the quadratic equation ( 131 ) for the complex frequency p. This 
equation can be used to determine p -for any given data. For the 
preceding data and for 


\s=0.7j ix = 1.0 


the natural frequency of oscillation was found to be I.287 cycles per 
second, with a logarithmic decrement of 1 . 370 . These results appeared 
quite satisfactory in regard to stability. 

It may be remarked that in the formulas of all the cases treated, 
the drag terms (induced plus parasite) had only a small influence on 
the stability characteristics of the rotor system. This shows the lack 
of necessity of determining the aerodynamic drag any more exactly than 
by the simplifying assumptions made in this ana lysis. 
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The differential equations of oscillation for any finite constant 
value of the speed ratio p e remain linear, as in hovering, hut they 
now have variable coefficients (periodic in i.e., in nt) instead 

of constant. In order to investigate the stability conditions in the 
transition from hovering to traveling, the speed ratio u was assumed 

to be a first-order small quantity. Solutions to the differential equa- 
tions could then be obtained by using the solutions in hovering as a 
first approximation and then making the corrections in accordance with 
the consistent procedure used in this report, that is, neglecting terms 
smaller than the second order. The correction in each case consisted 
in the addition of a particular integral of a non-hompgeneous linear 
differential equation with constant coefficients. These additions were, 
in all cases treated, found to be small in comparison with the corre- 
sponding solution in hovering. (A general development in powers of 
larger p e but < 1 might also prove convergent.) Cases A, B, and C 

of hovering were by this method treated for low-speed traveling, with 
the following results. 

In case A, the solution was found to be of the form 


0 = 0, 0 — 0^ + “ £2. 

where and ^ were the solutions for case A in hovering. In 
accordance with equations (99) end ( 108 ), 0-^ had the form 

, JL V 
% * Z,V 

n=»l 


where F was an arbitrary constant. The expression for 0 O was 
n *= 

then found to have the form 


n=l L- 


where 


A n and 


A * 


are constants which depend on F , 


given in general 


terms by equations (143a) and (143b) and for the numerical example by 
equations (144). The physical significance of these results is that 
four new natural frequencies appear, obtained by adding and subtracting 
fi/2ir to and from each of the two natural frequencies yj 2rj Q — and 

(l + rig) in hovering. It can be seen that three of the new frequencies 
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will not be significantly different from the havering frequencies, but 
that one new natural frequency appears which is approximately double 
the rotational frequency n /2it. The expression for P 2 also shows 

that the logarithmic decrements remain unchanged in the transition from, 
hovering to traveling. Therefore the solution which indicated no damping 


1 — " ft 

in hovering corresponding to the frequency \j2r\ e — still indicates a 
danger of instability in low-speed traveling. 


form 


In case B, the constraint condition was again assumed to be of the 


T «= *P 

As in the previous case, the solution for p was of the form 


3 - 3- + U~3„ 

1 e 2 

where P-j_ is the solution for case B in hovering. From equations ( 99 ) 
and (122) Is of the form 

- 2 pjif 

h-L !*»»“ 

n?=l 


The angle P 2 was then found to be of the form 


* 2 - 


= 5 


v 


(V 1 )* 

n 


where and A^’ are functions of F n given by equations ( 153 a) 

and (153b) in general terms. For the numerical example, A^F^ and 

shows 


n 


for 


A n */F n are given by equations (154). The expression for P 2 
that two new natural frequencies are added to the frequency O .775 — 

_ 2 rt 

Q ft 

hovering. These frequencies are 1.775 — and 0.225 — , neither of which 

2it 2 jt 

appears to present any particular danger of resonance. As in case A, the 
damping remains the same as in hovering; that is, the logarithmic decrement 
remains 1.282. This indicates satisfactory stability in the transition to 
traveling. 

In case C, the two kinematic constraints were also ass ume d as in 
hovering, namely. 
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(3 = X0, £ = (iS 

Proceeding as in the two previous oases, the solution for S’ was assumed 
in the form 


0 afl 1 + ^e s 2 

where 0 ^ was the solution for 


case C in hovering, and was given by 


Ea 

n=l 


P n t 


The constant B n is arbitrary and p n is given "by equation (131) in 

general and by equation (134) for the numerical example calculated. 
Then was seen to have the form 


^2 

n=l 


V 


(y 1 )* + A 


where and A n * are given as functions of B n by equation ( 165 a) 

and ( 165 b) generally, and by equation (l 66 ) for the numerical example. 

The expression for 8 ^ again shows that two new natural frequencies 

appear in the transition from hovering to traveling and that these 

can be obtained by adding (algebraically) ± ~ to the hovering natural 

frequency. For the numerical example treated (X « 0 # 7 , \x » 1) , the 

hovering frequency was 0.403 and the two new frequencies in low— 

2it 

speed traveling are therefore 1.403 -&• and 0.597 cycles per 

2jt 

second. The logarithmic decrement, as in cases A and B, remains the 
same as In hovering, namely 1.370 for case C. Thus the system appears in 
this case to remain quite stable in the transition from, hovering to 
traveling. 

It will be observed that the effect on the hovering oscillations 
of the transition to traveling is essentially the same for all the 
cases treated. 

In regard to the reaction of the blade system on the fuselage and 
to the strength of the hinge structure, the question of external and 
Intern al forces arising from the effect of external or internal constraints 
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is of interest. The external constraint moment caused hy the fixation 
of the pitch angle was given by the static and dynamic equations of the 
9 component, whereas the internal constraint moments act ing on the 
linkages "between the hinge axes were given hy the Lagrange multipliers 
of the derivatives of the kinematic conditions. Simple preadjustments 
were indicated between the angles G Q , and £ Q , by which the 

steady-state internal moments can be eliminated. 

In this paper the problems of vertical climbing, inclined travel 
direction, large speed ratios, disturbing external forces, and elastic 
vibrations have not yet been discussed. These problems can be solved, 
however, by the basic geometric, static, and dynamic procedure presented 
and by the introduction of the appropriate inflow velocities and inertia 
forces. 

Polytechnic Institute of Brooklyn 
Brooklyn, N. Y., May 5 , 19^6 
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APPENDIX A 


VELOCITY COMPONENTS IN PLANE OF CROSS SECTION OF BLADE 


By definition of the blade angle or pitch angle 9 , the line 

of zero lift makes an angle of 6 with the line h, where h is perpen- 
dicular to the z— aais and to the blade center line r. Let l , m, and n 
denote the direction cosines of any line with respect to the x— , y— , and 
z-axes, respectively. Then, from figure 1, the direction cosines of r 
are: 



Moreover, the direction cosines of line h, which lies in the xy— plane 
and is perpendicular to the projection of r on that plane, are: 


2^ = — sin £ 
o — cos £ 

“h = 0 

The direction cosines of c^ can now be determined as follows: Since 

c 1 is perpendicular to r, it follows from equation (Al) that 

l c ^ cos 0 cos £ — cos 0 sin £ + sin 0 = 0 (A3) 

Also, since the angle between c ^ and h is 0, it follows from equa- 
tion (A2) that 



cos 0 


2_ sin £ — m_ cos £ 
C 1 ‘ b c i 


(A4) 
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Moreover, 



+ n a 
G 1 


2 


= 1 


(A5) 


Equations (A 3 ), (A^), and (A 5 ) can. "be solved for the three unknowns l , m 

C 

an< -*' n 0n • ®' or Q sample, m and n can both be put In terms of 1 by 
x 1 1 c 2_ 

means of equations (Ah) and (A3). Substitution into equation (A5) then gives 
a quadratic in I c ^. The results are: 


= — cos 9 sin £ — sin p sin 9 cos 

m^ = — cos £ cos 9 + sin p sin 9 sin 

n « T cos 3 sin 9 
°1 + 


i 



(A 6) 


The lower alternative signs in parentheses, which really mean replacing 9 
by -0, must be rejected, because, as may be verified later (cf . equa- 
tion (10a) for L * ) , they would give negative lift in hovering. To 
« 

second-order small quantities, assuming in addition to relations (h) that 9 
is first order small, that is 


9 « 1 


equation (A 6) can be written as: 


Vl - - 5 + P0 


n - - 0 

C 1 


J 


(at) 


The oaarponents Y c ^ and Y n in the direction of c^_ and in the direc- 
tion of n normal to c^ and r will be: 
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= v Cl 


+ V. 




V n 
z c. 


Vn 


+ V. 


y 11 ^ 


V * n n 



(A 8 ) 


All quantities in equation (A 8 ) have already "been obtained except the 
direction cosines of n. These can be obtained by considering n as 
the vector product of unit vectors Cj^/fe i and v/r, inasmuch as n 

is perpendicular to both c^- and r. Hence 


= tn n — m n 
n c-j_ r r c-^ 


m. = n l - n I . 
n c i r r Ci l 


(A9) 




From equations (Al) , (A 7 ), and (A 9 ), the direction cosines of n are 
found to be: 


Z n » - (3 - £0 


m n “ — 9 + ^ > 



(A10) 


By putting equations ( 5 ), (A 7 ), and (A10) into equation (A 8 ), the expressions 
for V and V n given by equations ( 8 a) and ( 8 b) in the text are 

obtained. 
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APPENDIX B 

DETERMINATION OF LIFT COMPONENTS 


If it is assumed that the circulation vector F lies along_the_ 
direction r of the centerline of the "blade, the vector product r x Y 
(equation (6)) can "be written as follows: 

E * - - T a) + J - T =a) * - Vrjl < B1 > 


Therefore, from equations (5) end (Al) , the lift components per unit 
length will "be: 



By substituting the expression (9b) for the magnitude of T, equations (10a) 
of the text are obtained. 
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APPENDIX C 

EXPRESSIONS FOR HINGE MOMENTS 
The Moment M^ 


The first term in equation (22b) for the moment M^ can he expressed 
by the use of equation (8a) for V c ^ and equation (18) for c. Thus, 

if* V °V ■ f A3 °* 2 £ ls fe) m 


where 

V c \ 2 V 

~«M = s 2 + 2s|d e sin If + q e + n e C cos 
+ i n e 2 (l - cos 2i|r + 25 sin 2f) + 2^ e Ti e sin if 
to first-order small terms. 

Putting equation (C2) into equation (Cl) and integrating gives 

1 * a ±J=JL- 1 _ i] = ^—1 

3 


(C2) 


I. 


ds — f— = 1 + + + (1 + s*)^ sin * + r\ e + ii£ cos +) 


S i 1 - \& B 


+ hi 2 (1 — cos 2 1 + 2t sin 2if) + 2n n sin if - 

2 “ © 0 

” \ C 1 + s i + s i 2 )(^e sin * + ^e + ^ 008 ♦) 


(l + Sl )(l + Si 2 ) 


1 + s. 




"(1 — cos 2+ + 2i sin 2V) + U<i 0 T| e sin * 


•1 


= ~ [1 + s i + ^n e + 3d e 2 (l - Sj) + Ihx 0 (i + s 4 + 3^ Sin * 
+ lni e 5 cos t + & 0 ^ sin 2jr -'3H 0 2 (l - s^ cos aTj (C3) 
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when higher powers of s^ than the first and products of s^ and smaller 
quantities are neglected. 

From, equation (10b), the second term in equation (22b), to first- 
order small quantities, becomes: 



dr r 2 c 2 (l + (J. 


sin ^)Q + [i(3 cos y + 9 sin t — 7)J 


(CM 


where 


f 



Substituting equation (l8) for c, equation (C4), in the dimensionless 
variable s^= becomes: 

f f c dr =» fjtpn^-^c. 2 P ds 1 ' ~ ( s + sin [ss + u (P cob i(r 

i J 3 ! 1 - 13 ! L 

+ 0 sin * - 7 )j « f 7tp^ 3 Ci 2 ^(l + s*) - u e g(l + s ± ) + d e |(i + s i) e 

- - Si )] sin y + £(l + cos t + “* s j) !3 sln 

+ -~-(l “ s^ 9(1 — cos 2*)|- 


2 * 


(C5) 


Addition of equations (C3) and (C5), together with the third term 
of equation (22b) yields the expression for given by equation (30a) 

in the text. 


The Moment 

By use of equations (10b) and (l8), the first term in the expression 
(equation (23)) for can be written as follows, to first-order small 

quantities: 
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9 ^ 




L *rdr = — rtpf 


r ds — — ? ■ [s 3 e + s 2 u.(3 cos i|r + 20 sin * 

\k \/^7 L 


+ su e 2 (]3 cos ^ + 0 sin T|f - 7 ) sin 


( 06 ) 


The integral in equation (C6) can he evaluated by means of the change 
of variables 


\J 1 — s « u, ds = — 2u du, s = 1 — u c 


and the result is given hy equation (31a) of the text. 


The Moment M ^ 

The moment M ^ as given hy equation (28a) will he second order 

small. Hence it is necessary here to use terms of the second order. 
From equations ( 10 a) and (lh). 



where A — £ cos tJt + e sin 'if — y. Thus, hy using equation (18), the 
first term of equation (28a) for M zl becomes: 



where M and IT are the coefficients of u and (i 2 , respectively, in 
equation (C7) . The integral of equation (C8) can he evaluated in the 
same way as that of equation (C6), and the result is given hy equa- 
tion (32a) in the text. 
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APPENDIX D 

INCREMENTS IN LIFT COMPONENTS DOE TO DAMPING 

From equation (9a), tlie total circulation (i.e., including the 
velocities |3 and £ ) will "be 


r “ KO (?n + A 7 n) 


(Dl) 


where V is given "by equation ( 8 b), and AY is the increment in Y n 

v* * a 

due to p ertfi ^ . This increment can be obtained from the relation 
(cf. equations (a8 )) 


Y = l_ AT_ + m AT + n AT (D2) 

n n x n y n z 

where, for example, AT « T . + V .. Thus, by using equations ( 38 ), (39), 

x xP x£ 

and (A10), the expression for AY n becomes: 


0 


which to second-order small terms reduces to 

AT • « 

-£ = -« i - £■ o>3) 

Gr 0 n 

The total circulation will therefore be 

r= Jtcft* |^ 0 (l + n) + h£(3 + £ 0 ) cos + + (® - 35) sin Y—t| — ® 

From the Eutta-Jouhowski relation (equation (Bl)), the direction 
cosines (equations (Al)) of the line r, and the velocity increments 
(equations ( 38 ) and (39)), it is seen that the expressions (equations (B2)) 
for the lift components per unit length must be changed by the following 
increments: 



96 


HACA TN No. l^O 


£L • 
X 

prnr 


L 3 _ £ £ = AK 
n « x 


pFS* 



L * 


z 


pTfZr 




(I>5) 


The total lift components, for example the x-components, may then he 
written as follows: 


V “ pflr(r + Ar)^ + 


0 56) 


where r is given hy equation (9b), AT (see equation (D3)) hy 


and k hy 



(D7) 


to he taken from equation (B2) . By neglecting cross products of damping 
terms, the additional lift components can, according to equation (D 6), 
he written in the form 

(d8) 

Thus, hy the use of equations (B2), (l>7), (9b), and (D5), equation (D8) 
leads to the results given hy equations (^la), (Ulh), and (ilc) in the 
text. 
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APPENDIX 1 

DAMPXNG--MOMEMP INCREMENT ABOOT BLADE AXIS 

The following relations follow readily from figure 6 , where the 
velocities (incl uding 6s in the vector diagram) shown are those of the 
relative wind: 



where V can he taken from equation (8a), and 
C 1 
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The integration gives 

M “ ~ li P®*® 2 ®! 3 “ + s i)(f + s i + 7n e + 7n e (sin if + £ cos +)| 
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APPENDIX F 


DAMPING-MOMENT INCREMENTS ABOUT BINGE AXES 


Increment in M 


xl 


In order to evaluate the first term of equation (46), the increment 
AY Ci must first he determined. This can he easily obtained from 

(cf. equation (A8)) 


AY = 7 AY„ + m. AY + n AY 
cq °i x c i y c i 2 


(FI) 


The direction cosines !„ , m , > and n are given hy equation (A7) , 

°1 C 1 C 1 

and the velocities A Y x , and so forth are given hy the addition of equa- 
tions ( 38 ) and (39). Thus, it is found that up to first— order small 
quantities 


AY„ 


-fir L 


Q 


(F2) 


(Since C M is already first order small, it suffices to obtain AV 

ac C 1 

to only first orders in order to determine the first term of equation (46) 
up to second— order small quantities, ) 

1 

By using equation ( 8 a) for Y .and equation ( 18 ) for c it is 

°1 

seen that to first— order small terms 
f dr c 2 Y Cl AY c _^ = — ^ P ds 1 — S - (s 2 + sp.^ sin 

J r ± 11 J a 1 - b ± 


- - 0^ \ 2 1 ( x * s i) C 1 - % aln ♦) < F3) 


to first powers of s^. 

From equations ( 18 ) and (4lc) the second term of equation (46) is 
seen to he: 
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Jr t v ° dT ’~’ v i : dr (1 + T| + n£ cos + u Bin if) 

+ A ^20 + |i((3 cos * + 9 sin if - r)J 

/ S ' as (rr^<| 8a + « AS + « B3 


p “3 p 

= - Jtpfi xl°o^ 


+ 2 i 


0S 


(Fit) 


ft 


where, temporarily, the following abbreviations for terms, not containing 
the integration variable s (or r), have been introduced: 


A s ti + u (sin t + 5 cos if) 

'G © 


D s n (p cos i|r + 0 sin if - 7) 
© 


(F5) 


The integral in equation (Fit) can be easily evaluated, with the result : 

,(£ cos t + sin ^ 


r AL z *c <3r a - *pfl% 3 Cl 2 (l + °i) ’ I JT O’ + ^e + 2u e< 

L 

+ ft J [f + 008 + + 9 sin t - 7) j 


(F6) 


The third term of equation (46), with the substitution of equa 
tion (8a) for F Cl and equation ( 18 ) for c is: 

a T 0i c3 = | *V (s * A)|n i-J /2 


ds (F7) 


where A has the same meaning as in equation (F5) • 

The integral in equation (F7) can be evaluated by means of the 
change of variables 
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1CBL 


1 — s = u, els 


du. 


S a 1 — U, 



(18) 


Thus, . 

6 f V c 3 dr » £ n% 2 c 3 f 1 (1 + A - u) ■-- - •■ ■ du 

Jr°l a ij n «3/a 


" | fl2R2 °i 3 jj“ ( x + A ^( 1 “ s i) 

a — fl^R 2 ^ 3 — (l + e ^[2 + s 

n 1 35 V iV{_ 1 

+ 7^1 — 2s ^)^ e ( sin t + £ cos 



(F9) 


By putting equations (F 3 ) , (F 6 ), and (F 9 ) into equation (46), the result far 
the damping moment as given "by equation (4-9) of the text is obtained. 


Damping Moment 

Putting equation (i^lc) into equation (4-7) shows the expression for the 
moment to he: 




o U 1 
r dr = JtpflR 


r 1 Y^Y 2 

J b i v 1 - «i / 


,r. n ) 


+ s 3 ^i + 2 £■ ej 


(F10) 


where A and D are given by equations (F5) ■ The integrals in equation (F10) 
can be evaluated by means of the change of variables (equation (F 8 )), with 
the results: 
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^ “ s l. 


+ §0 

105 \ 2 V 


Z 11 d8 /in±_y s = 3 - 22_ Z + ie ') 
t/ 8 ! V; - \) 315 V 2 V 


The use in equation. (F10) of these values leads to the expression (aqua — 
tion (50)) for the damping moment AM^ given in the text. 


Damping Moment AM 


zl 


According to equation (48a), 

B 


AM. 


zl 


/l E p E 

AL *r dr + / AD *r dr 

r l U r i y 


(Fll) 


From the expressions for AL * (equation (4lb) and for c (equation (18)) 
the first term of equation (Fll) is seen to he 

P ALy*r dr «= — 


/ mi (^y ie + s2|j ° (ap °°° * 

+ 0 sin t 


With the preceding values of the integrals, equation (F12) becomes: 
Ayr dr . - ♦ |ej) | 


(F12) 


32 15 

— “ 9 + W_(2p cos + + 0 sin + 

315 105 e 


“ 2 7> 


(F13) 


From the expression (equation (44)) for AD * , the second term of 
equation (Fll) is seen to be 


/. 


E 


AD *r dr = 


R 


1 t M + 9 

r . ^ i°p + 5 " ) (1 + H sin 1f)ncpn 2 r^ 

5 1 + H sin f J ft 


(F14 
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where D is given by equation (F 5 ) . Substituting equation ( 18 ) for c, 
integrals s imilar to those in equation (F12) appear, and the result is: 


£Dy*r dr = h p S?r\(2. + ^i) £ I ^ 


l 2 a + - 6 
< P 5 


+ 3li la, sin + i(p 00 s t + 8 sin t — 7 ) ► (F15) 
e|_P 5 J 


By addin g equations (F13) and (F15), the expression (equation (51)) for the 
damping moment £M Z -^ is obtained. 
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APPENDIX G 

PHYSICAL CONSTANTS OF BLADE SYSTEM 
General Assumptions 

The following is a summary of the data assumed in this report for the 
purpose of numerical calculation: 

r < 

a = -i. a 0.2 
i R 

tv n 

— b Constant along span 1 of blade * — 
c 8 

n b 4 blades 

£ « 0.0025 

a ■ 0.020 
P 



W = 4000 lb 

R = 25 ft 

Cl = 20 radians/sec 



Moments of Inertia 

Because the pitch angle, with the foregoing data, will "be ijuite small 
( 0 C = O.O 306 , equation ( 103 )) it will be sufficiently accurate, for the 
purposes of calculating the moments of inertia, to neglect the rotation 
of axes due to the angle 0 . 

In general, for a blade, 

I max) = f l°\' x 2 ( E X min) = W 0 (01) 

where and f are constants, which for a solid Clark Y section have 

the values 


= 0.0418, f 2 = 0.0454 
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From, equation (Gl) , 


hJjLfi .f 

II f x W 0.0413 \SJ 


0.017 


(G?) 


Neglect, for simplicity, of I p in comparison with 1 ^ will therefore 

he justified. Similarly, moreover, in accordance with the definitions in 
equations ( 6 l), the moment of inertia I may he neglected in comparison 

with I__. It follows then, because for an airfoil section such as 

Clark Y, I, and I will he practically the same, that 
J. zs 

(fe S + z s 2 ) “■ * X 1 + h * \ - J S * 

Therefore, with the notations in equations ( 6 l), 


I 


r 




(G3) 


The numerical value of I can h? readily determined hy using equation (Gl) 
Thus 



The value of hy definition, will he 

%-r 


Ar 2 dr 


(G5) 


where the cross-sectional area A may he given hy 



m h, | h| 
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and f^ is a o oust ant, which for a solid Clark Y section has the 


Hence 


f* 3 = O.725 




E V r 2 dr f R c 2 r 2 dr 

% 8 Jr ± 


‘3 2 /' R ^ 2 “ R ^ _ f 3 2 C 1 + S i) 

r °i 


f r 2 — — - dr - 12 . e3 c 2 

(K c- 1 -B) te " r 01 12 


i!i B 4 (j Ifl) 

6 8 1 12 


From equations (Gk) and (G 7 ), 

,3 


a 24f. 




2k x 0.0kl8 x 0.6 x [ y ® 0.00386 


Therefore I may "be considered negligible in comparison with 1^. 
From equation (G 7 ) , it follows that 


R 


‘<=1 - 7T 0- - = 0 % - ZZZ C 1 - * 0 % 


0.725 


796 (1 - B>f. 


■H 


The value of S is, by definition. 


1: 


Ar dr 


value 


(07) 


( 08 ) 


(09) 


(010) 
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Putting equation (G 6 ) into equation (G10) shows the value of § to he: 


?? 1 3 p 3 

5 z r 1 


(1 + 8i ) 


(Gil) 


Comparison of equations (G 7 ) and (Gil) shows that 

S = 2^ 

R 


(G12) 


Thus all the relations in equations (77) of the text have been derived. 
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APPENDIX H 


EXPLICIT DERIVATION OF OSCILLATIONS IN CASE C 


By putting equation (126) into equations (97 a ).» (97b) * and (97 °) > 
the expressions for the P*s, in accordance with the definitions 
(equation (128)), are found to he the following: 


where 
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Equation (130b) therefore becomes: 


p 2 ! + Hp — J + X 


E_ I R 1 * 0 2id 

n 2 12 E V 


C 1 + 


- x8 [»% + 9 ^ ^ ( x ♦ !"•) *^ + a 0 
! (p% * 2n e %) + ** |R? 2 8 c i - (1 * |«i) H5 (| e o " a P )» - 0 


or, rearranging in' power s of p and. dividing through by I. gives 


p 2 ( — X 2 - n 2 + i- + 2 X 41 0 


1 + s, + 2nj o. 


x2\ +p |1. + ^(m -+ 3^ E / o\ 

V K b 315 : 


f II - X 2 




315 




- Xn 


s 32 

«% 315 


(j + K)$ v«p) 


L_xa(i + a g ♦ * y ♦ |*i) “ “V* a _ “ 


Equation (H2) is eqxd.vaJ.ent to equation (131) of the text, when 

has been tahen as negligible in comparison with 


(x 2 + H 2 ). 

Equation (H2) is most easily handled by using numerical values 
immediately. With the numerical data consistently assumed in this report 
(see appendix G-) and with the use of equations (77), the following values 
are obtained, in accordance with the abbreviations in equation ( 131 ) : 
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E 32 

315 

a h 


(l +■ is ( ^0 — a^X — a — 

\ 2 y^5 ° ■ p; p 5 c 


^4 

= ir § 1 32_ /l.l X 0.00U5X - 0.020 - - X 0.0301) 

0 T_ 315 V 5 / 


I H 315 


32 


m 3.14 x 0.0025 x 796 x 0.8 x — (0.00495X - 0.032) 

315 


0.00252X - O.OI63 


R^c , 


«- . - t % (till 111 (| x 0.15 x 0 

in ° 12 \ V 6 


.0306 - °^£) 

6« J 


- 3.14 x ° - , - 99 - 2 ^ x 796 x 0.8 x O.OO165 - - 0.00069 


5- = 3.14 x 0.0025 x Z 2 6 ..- *-P - ^ x i— x — - x 2.55 - 0.00168 
X 12 36 35 


H 


~ * 3.14 x O.OO25 x ? - 2 6 - — l 8 - x i x 0.15 - O.OO386 = 0.00656 

\ 12 6 


E 32 


(l + * 3-14 x 0.0025 x 796 x 0.8 x x I.075 « 0. 


548 


^3!5 

E fl + s. + 2n } c. 10. 

—5—- -i = — f — » 3.14 x 0.0025 X 796 x 0.8 x — !-■ x 0.15 X i - 0 

£1TT_ 12 R 12 6 

H 


Substitution of the foregoing numerical values into equation (H2), the 
changing of all signs there, and neglect of the quantity 

*1 J \ / p 9 \ 

— + 21ji0 c — I In comparison with — l \ + u J leads to the following 

* v 

equation: 


.0136 
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■/ 

p 2 (x 2 + n 2 ) + p( 0.00252x4-1 - 0 . 01631 a + 0.0006911 - 0.00168 

w 

+ 0 . 5981 2 - O.OI 36 X) 

+ ( 0.00656 + 1 . 1 A 2 + O.lii 2 - 0 . 5601 ) a 0 (H 3 ) 

This equation is equivalent to equation (132) of - the terfc. 


* 




Figure 1.- Blade position. 



2.~ Components of traveling velocity v. 






(b) Plan, view. 


Figure' 3.- Centrifugal force element on 


blade and axes x-^, y, , of moments 
about hinges. 



Figure 5.- Lever [r(f; - ^)J of centrifugal force 

*(£-£].) e p 

element — or (£ - CjJs £§ = £,. 
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